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ABSTRACT 

It has long been known that photoionization, whether by starlight or other sources, has diffi- 
culty in accounting for the observed spectra of the optical filaments that often surround central 
galaxies in large clusters. This paper builds on the first of this series in which we examined 
whether heating by energetic particles or dissipative MHD wave can account for the observa- 
tions. The first paper focused on the molecular regions which produce strong H2 and CO lines. 
Here we extend the calculations to include atomic and low-ionization regions. Two major im- 
provements to the previous calculations have been made. The model of the hydrogen atom, 
along with all elements of the H-like iso-electronic sequence, is now fully nZ-resolved. This 
allows us to predict the hydrogen emission-line spectrum including excitation by suprathermal 
secondary electrons and thermal electrons or nuclei. We show how the predicted H I spectrum 
differs from the pure recombination case. The second update is to the rates for H** - H2 inelas- 
tic collisions. We now use the values computed by Wrathmall et al. The rates are often much 
larger and allow the ro-vibrational H2 level populations to achieve a thermal distribution at 
substantially lower densities than previously thought. 

We calculate the chemistry, ionization, temperature, gas pressure, and emission-line spec- 
trum for a wide range of gas densities and collisional heating rates. We assume that the fil- 
aments are magnetically confined. The gas is free to move along field lines so that the gas 
pressure is equal to that of the surrounding hot gas. A mix of clouds, some being dense and 
cold and others hot and tenuous, can exist. The observed spectrum will be the integrated emis- 
sion from clouds with different densities and temperatures but the same pressure P/k — riT. 
We assume that the gas filling factor is given by a power law in density. The power-law index, 
the only free parameter in this theory, is set by matching the observed intensities of IR H2 
lines relative to optical H I lines. We conclude that the filaments are heated by ionizing parti- 
cles, either conducted in from surrounding regions or produced in situ by processes related to 
MHD waves. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: NGC 1275 - galaxies: 
clusters: individual: NGC 4696 - intergalactic medium - infrared: galaxies 



1 INTRODUCTION 

The central galaxies in large clusters are frequently surrounded 
by a system of filaments that emit strong molecular, atomic, and 
low-ionization emission lines. Understanding the origin of this 
line e mission has been a long-standing challenge ( Johnstone et al 
l2007h . The line luminosities are too great for the filaments to be 
powered by known sources of radiation such as the central AGN 
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or the diffuse emission from the surrounding hot gas. The evolu- 
tionary state of the gas is totally unknown. Two possibilities are 
that they form from surrounding hot gas or by ejection from the 
central galaxy in the cluster. Star formation may occur in the fila- 
ments although the emission-line spectra do not resemble those of 
Galactic H II regions near early-type stars. Questions concerning 
the origin, energy source, and evolutionary history are important, 
in part because of the large mass that may be in volved, as high as 
~ 4 X 10'° Me according to lSalome et all jzOOd) . 

Infrared spectra show H2 lines which are far stronger rela- 
tive to H I lines than those emitted by molecular gas near O stars. 
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Attempts at reproducing the spectra assuming starlight photoion- 
ization are largely unsuccessful, as reviewed by iJohnstone et all 
l l2007h . Hybrid models, in which starlight produces the optical 
emission while other energy sources produce the molecular emis- 
sion, appear necessary. The fact that optical and molecular emission 
luminosities trace one another dJaffe et alll2005[) would require that 
these independent energy sources have correlated luminosities. 

Since photoionization by O and B stars does not app ear able 
to account for the spectrum, Paper 1 jperland et ail l2008) consid- 
ered whether purely collisional heating sources can reproduce the 
observed spectrum. Magnetic fields occur in the environment and 
wave energy is likely to be associated with the field. Dissipation of 
this wave energy could h eat the gas. Cosmic rays are also present 
dSanders & Fabianll2007h . These would heat and ionize the gas and 
produce strong low-ionization emission. MHD waves may also ac- 
celerate low-energy cosmic rays within a filament. Similar particle 
and wave processes occur in stellar coronae where they are collec- 
tively referred to as non-radiative energy sources, a term we shall 
use in the remainder of this paper 

Paper 1 focused on H2 lines produced in molecula r regions 
that are well shielded from light. Ijohnstone et all i2007h detected 
H2 lines with a wide range of excitation potential and found a cor- 
relation between excitation and the derived population temperature. 
We found that non-radiative heating with a range of heating rates 
can reproduce the observed H2 spectrum. In this paper we con- 
centrate on atomic and low-ionization emission and develop the 
methodology needed to predict the spectrum of gas with a range of 
densities but a single gas pressure. We find that cosmic-ray heating 
produces a spectrum that is in general agreement with a wide range 
of observations. Purely thermalized energy injection cannot repro- 
duce the spectrum. This does not rule out MHD wave heating but 
does require that they produce or accelerate high-energy particles 
in addition to other forms of energy. The resulting model, while 
empirical, points the way for a physical model of the origin of the 
filament emission. 



2 SPECTRAL SIMULATIONS 
2.1 The basic model 

Starlight photoionization has long been known to have difficulty 
in reproducing observations of cluster filaments. Hybrid models, in 
which different energy sources produce the molecular and atomic 
emission are more successful but have problems accounting for 
why different energy sources would correlate with one another. 
Here we assimie that only non-radiative heating, by either cosmic 
rays or dissipative MHD waves is important. The entire spectrum 
is produced by these energy sources. 

For simplicity we assume that the gas is well shielded from 
significant sources of radiative heating. In reality light from the 
central galaxy or the surrounding hot gas will photoionize a thin 
skin on the surface of a cloud but will have little effect on the ma- 
jority of the cloud's core. The emission lines emitted by the ionized 
layer will be faint. In the calculations that follow we include the 
metagalactic radiation background, including the CMB, so that the 
continuum is fully defined from the gamma-ray through the radio. 
As described in Paper 1 this external continuum is extinguished 
by a cold neutral absorber with a column density of 10^' cm^^ to 
approximate the radiation field deep within the filaments. This con- 
tinuum is faint enough to have little effect on the predictions in this 
paper. 



Table 1. Assumed gas-phase abundances by log nucleon number density 
[cm^^] relative to hydrogen. 



He 


-1.022 


Li 


-10.268 


Be 


-20.000 


B 


-10.051 


C 


-3.523 


N 


-4.155 





-3.398 


F 


-20.000 


Ne 


-4.222 


Na 


-6.523 


Mg 


-5.523 


Al 


-6.699 


Si 


-5.398 


P 


-6.796 


S 


-5.000 


CI 


-7.000 


Ar 


-5.523 


K 


-7.959 


Ca 


-7.699 


Sc 


-20.000 


Ti 


-9.237 


V 


-10.000 


Cr 


-8.000 


Mn 


-7.638 


Fe 


-5.523 


Co 


-20.000 


Ni 


-7.000 


Cu 


-8.824 


Zn 


-7.6990 







In keeping with our assumption that the regions we model are 
well shielded, resonance lines such as the Lyman lines of H I or the 
Lyman- Werner electronic systems of H2 are assumed to be opti- 
cally thick. Because of this continuum fluorescent excitation of H" 
and H2 does not occur. With these assumptions the conditions in the 
cloud are mainly determined by the non-radiative heating sources 
which are the novel aspect of this work. 

For simplicity we assume that the chemical composition is 
the same as the Local ISM. The detailed gas-phase abundances are 
based on emission-line observations of the Orion star-forming en- 
vironment and are given in Table [T] A Galactic dust-to-gas ratio is 
assumed. Refractory elements are depleted from the gas phase in 
keeping with our assumption that dust is present. 

Grains have several effects on the gas. H2 forms by catalytic 
reactions on grain surfaces in dusty environments. Collisions be- 
tween gas and dust tend to bring them to the same temperature. 
This process either heats or cools the gas depending on whether the 
grain temperature is above or below the gas temperature. Molecules 
can condense as ices coating the grains if the dust becomes cold 
enough. Each of these processes is considered in detail in Cloudy, 
the spectral synthesis code we use here, but the underlying grain 
theory depends on knowing the grain material, its size distribution, 
and the dust to gas ratio. The grain temperature depends on the 
UV-IR radiation field within the core, which in turn depends on 
whether in situ star formation has occurred. Rather than introduce 
all of these as additional free parameters we sim ply adopt t he grain 
H2 catalysis rate measured in the Galactic ISM t Jurall 19751) . We do 
not consider grain-gas energy exchange and neglect condensations 
of molecules onto grain ices. Tests show that these assumptions 
mainly affect the details of the H" - H2 transition. One goal of our 
work is to develop a physical model that accounts for the spectral 
observations. A long-term goal is to use such a model to determine 
the grain properties from observations of line extinction and the 
infrared continuum. 

Substantial uncertainties are introduced by the need to assume 
a specific gas-phase composition and dust properties. It would be 
surprising if the gas and dust composition happened to match that 
of the local ISM, although it would also be surprising if it were 
greatly different. The molecular collision rates, described in Sec- 
tion l2.4l have their own substantial uncertainties, probably 0.3 dex 
or more. These considerations suggest that there is roughly a factor- 
of-two uncertainty in the results we present below. This is intended 
as a 'proof of concept' calculation aimed at identifying what phys- 
ical processes may power the observed emission. If successful, we 
can then invert the problem and determine the composition or evo- 
lutionary history from the spectrum. 
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2.2 Non-radiative energy sources 

Paper 1 considered two cases, heating by dissipative MHD wave 
energy, which we assumed to be deposited as thermal energy, and 
cosmic rays, which both heat and ionize the gas. We refer to these 
as the 'extra heating' and 'cosmic ray' cases below. The effects of 
these energy sources on the microphysics are fundamentally differ- 
ent. 

Supersonic line widths, often thought to be due to MHD waves 

associated with the magnetic field, are observed in the ISM of our 

I — 1' 1 

galaxy ( IDvson & W illiams 1997). The field and wave kinetic en- 
ergy are often in rough energy equipartition ( Mvers & Goodman 
Il998t) although the details remain uncertain. iHeiles & Crutchei 
yOOa) review the extensive numerical simulations of MHD waves. 
Waves can be damped by processes such as ion-neutral drift, 
which convert wave energy into other forms of kinetic energy, al- 
though the details are uncertain and the process may be unstable 
jTvtarenko etaill2002h . In this exploratory work we simply want to 
quantify the effects of such heating on otherwise well-shielded gas. 
We adopt a heating rate that is proportional to density. 



actions with low-density gas ar e describ ed in Spitzer & Tomasko 



1968 ), iFerland & Mushotzkvl 1 1984 ). lOsterbrock & Ferlanc 



, |Fm 

2006h. IXu & Mc Cravl h99lK Ixu & McCravl ( Il99lh 



// = /fo[«(H)/no] [ergcm-^s-'] 



(1) 



This form, modified from that used in Paper 1, includes the ratio of 
the hydrogen density n(H) to a scale density no which we take as 
1 cm^-' . This density dependence causes the wave-heating rate to 
go to zero as the gas density goes to zero. 

We parameterize the extra heating rate as the leading coeffi- 
cient in the heating rate H(). We assume that the heating simply 
adds to the thermal energy of the gas so that the velocity distribu- 
tion remains a Maxwellian with a well-defined kinetic temperature. 
With these assumptions the only collisional processes which occur 
are those which are energetically possible at the local gas kinetic 
temperature. This is a major distinction between the extra-heating 
case and the cosmic -ray case described next. 

The second case we consider is energy deposition by ioniz- 
ing particles. These particles could be related to the high-energy 
particles which are known to exist in the hot gas surrounding the 
filaments dSanders & Fabianll2007l) , or could be caused by MHD - 
related phenomena like magnetic reconnection jLazaria nl l2005h . 
Whatever their fundamental source we will refer to this as the 
cosmic-ray case for simplicity. As in Paper 1 we specify the ion- 
izing particle density in terms of the equivalent cosmic-ray den- 
sity relative to the Galactic background. We adopt t he background 
cosm ic ray H? dissociati o n rat e of 3 x 10^'^ s^' jWilliams et all 
1 1998b . ISanders & FabianI ( 12007 ) find an electron energy density 
that is roughly 10^ times the galactic background in inner regions 
of the Perseus cluster. This value guided our choice of the range of 
cosmic-ray densities shown in the calculations which follow. 

There is some evidence that this Galactic backgro und cosmic 
ray H 2 dissociation rate may be substantially too low. I Shaw et all 
( I2OO8I) found a cosmic-ray ionization rate forty times higher 
along the sight line to Zeta Per from detailed modeling while 
llndriolo et al ( |2007|) derived a value ten times higher from the 
chemistry of H3" along fourteen different sight lines. If these newer, 
substantially higher, Galactic background rates are accepted as typ- 
ical then the ratio of the cosmic ray rate to the Galactic background 
given below would be reduced by ab out an order of m agnitude. This 
is clearly an area of active research (Dalg arndl 19961) . We adopt the 
Galactic cosmic-ray background quoted in Paper 1 for consistency 
with that work. We will express the particle ionization rate in terms 
of the Galactic cosmic ray background rate. 

Cosmic rays both heat and ionize the emitting gas. Their inter- 



iDalgar no etal (1999)rTi neetall ^1997^ and many more pa- 
pers. [Abel et al ( 2005 ) describe our implementation of this physics 
in the current version of Cloudy. Briefly, if the gas is ionized (the 
electron fraction iie/nu > 0.9) cosmic rays give most of their 
energy to free electrons which are then thermalized by elastic 
collisions with other electrons. In this highly-ionized limit cosmic 
rays mainly heat the gas. In neutral gas, with low electron fraction, 
some of the cosmic-ray energy goes into heating the gas but much 
goes into the creation of a population of suprathermal secondary 
electrons. These secondaries cause direct ionization as well as 
excitation of UV resonance lines. These excitation and ionization 
rates depend on the cosmic-ray density and the electron fraction 
but not on the kinetic temperature of the thermal gas. 

These two cases behave quite differently in the cold molec- 
ular limit and in the nature of the transition from the molecular 
to atomic and ionized states. In a cosmic-ray energized gas the 
suprathermal secondary electrons will cause ionization and excite 
resonance lines even when the gas is quite cold. As the cosmic-ray 
density is increased the transition from the fully molecular limit 
to the atomic or ionized states is gradual since the ionization rate, 
in the low electron-fraction limit, is proportional to the cosmic-ray 
density and has little dependence on temperature. Thus a very cold 
cosmic-ray energized gas will have a significant level of ioniza- 
tion, dissociation, and excitation of the H I, He I, and H2 resonance 
transitions. The UV lines, being optically thick in a well-shielded 
medium, undergo multiple scatterings with most lines being ab- 
sorbed by dust or gas. Some emission in infrared and optical sub- 
ordinate lines, and in ro-vibrational transitions in the ground elec- 
tronic state of H2, occurs. 

In the extra-heating case the thermal gas has a single kinetic 
temperature. As the extra-heating rate increases the temperature 
rises but the gas remains fully molecular until the kinetic temper- 
ature rises to the point where dissociation is energetically possi- 
ble. The transition from molecular to atomic phase occurs abruptly 
when the gas kinetic temperature approaches the H2 dissociation 
energy. The gas will be in one of three distinct phases, with es- 
sentially all H in the form of H2, H*', or H+. As the heating rate 
and kinetic temperature increase the gas will go from one form to 
another in abrupt transitions which occur where the temperature 
reaches the appropriate value. 

These two models mainly differ in how energy is deposited in 
molecular or atomic gas. The energy of a cosmic ray will eventually 
appear as heat, internal excitation, or as an ionizing particle, while 
in the wave-heating case the energy deposition is purely as heat. 
The cosmic-ray case does not exclude MHD effects as the funda- 
mental energy source since MHD waves can produce low-energy 
cosmic rays within a filament ( lLazariadl2005l) . 

2.3 A fully /-resolved H-like iso-sequence 

The calculations presented here use version C08.00 of the spec - 
tral synthesis code Cloudy, last described bv lFerlande"tail ( ll998h . 
There have been several improvements to the simulations since Pa- 
per 1. These are described next. 

The hydrogen recombination spectrum is special because it 
can be predicted with great precision. Two of the most important 
rates, radiative recombination and transition probabilities between 
bound levels, are known to a precision of several percent. If the H I 
lines form by recombination in an ionized gas that is optically thick 



© RAS. MNRAS 000.mi27l 



4 



G.J. Ferland et al. 



in the Lyman conti nuum then the emission -line spectrum should be 
close to Case B of iBaker & Menzell ( Il938ll . Collisional excitation 
of excited levels from the ground state of H" are a complication 
since their contribution to the observed H I spectrum depends on 
density, temperature, and ionization in ways that are fundamentally 
different from the recombination Case B. 

The computed H I emission-line spectrum presented here is 
the result of a full calculation of the formation of the lines, in- 
cluding radiative and three-body recombination, collisional ioniza- 
tion, collisional coupling between nl terms, and spontaneous and 
induced transitions between levels. For much of its history Cloudy 
has used the compact model of the H-like iso sequence described 
bylF erguson & Ferland ( 1997) with the rate formalism described by 
iFerla nd & Rees ( 1988). That model resolved the 2s and 2p levels 
but assumed that /-terms within higher-« configurations were popu- 
lated according to their statistical weights. This is formally correct 
only at high densities. Various approximations were introduced to 
allow the model to work well in the low-density limit. 

Our treatment of the entire He-like iso-electronic sequence 
was recently enhanced to resolve any number of nl levels, resulting 
in a far better representation of the physi cs of the re s ulting emis- 
sion. Application s to He I are discussed bv |Porter"et ail( l2005l) while 
IPorter & FerlmS bOOTh describe ions of the He-like sequence. 
This same methodology has now been applied to the H-like iso 
sequence. Levels can now be fully /-resolved and any number can 
be considered. The m ethods used to cal culate rates for /-changing 
collisions are given in FPorter et all ( |2005|) . 

A very large number of le vels, and resulting high precision, 
were used in IPorter et all feOOSh . The was possible because of the 
overall simplicity of that study. The temperature and electron den- 
sity were specified so a single calculation of the level populations 
was sufficient to obtain the He I spectrum. The calculations pre- 
sented below self-consistently solve for the temperature, ionization, 
and chemistry, which requires that the model atom be reevaluated 
a very large number of times for each point in the cloud. An in- 
termediate model must be adopted with a small enough number of 
levels to make the problem solvable with today's computers but yet 
expandable to have more levels and become more accurate when 
more power becomes available or high precision is needed. A hy- 
brid approach was chosen. The lowest n' configurations are nl re- 
solved. They are supplemented with another n " configurations, re- 
ferred to as collapsed levels, which do not resolve the nl terms. 
The number of resolved and collapsed levels can be adjusted with a 
tradeoff between accuracy and speed. In the following calculations 
levels ;j ^ 15 are nl resolved with another ten collapsed levels rep- 
resenting higher states. This is sufficient to achieve a convergence 
accuracy of better than 2% for the H I lines we predict. Larger er- 
rors are introduced by uncertainties in the collisional rates. These 
are difficult to quantify but probably produce errors of ~ 20% in 
the line emissivities and ~ 5% in relative intensities of H I lines. 

Resolving the nl terms makes it possible to predict the de- 
tailed effects of thermal and cosmic-ray excitation of optical and 
infrared H I lines because the resulting emission spectrum is sen- 
sitive to the precise nl populations. High-energy particles will 
mainly excite np levels from the ground l.s term because most 
atoms are in the ground state and electric-dipole allowed transi- 
tions hav£j^ge££ollisional_excitation cross sections at high ener- 
gies jSpitzer & Tomaskoll 19681) . After excitation, the np level will 
decay to a variety of lower j or d terms because of selection rules 
for optically-allowed transitions. The H I Lyman lines are assumed 
to be optically thick so that np—t Is transitions scatter often enough 
to be degraded into Balmer lines plus La. The resulting La pho- 



tons will largely be absorbed by grains rather than emerge from the 
cloud. 

State-specific cosmic-ray excitation rates to np from ground 
are derived from the Bom approximation. This is appr opriate due to 
the high energy of a typical seco ndary electron |Spitzer & Tomasko 
1968). The theory described in lAbel et all ilOO^ and Ishaw et al 
(l2008) is used to derive a secondary excitation rate of La, 'fif_^2p' 
The secondary excitation rate of any permitted line, is given in 
terms of the La rate by 

where £ is the energy of the transition and gf is its oscillator 
strength jShemanskvetall 19851) . 

Quantal calculations of collisional excitation rates for exci- 
tation of nl terms by thermal electrons are used jAnderson et all 
|200C|) . There is no favored final / term, another distinction between 
the cosmic ray and extra heating case. Details of the resulting opti- 
cal and infrared H I emission will be discussed below. 

2.4 Revised H - H2 collision rates 

The rates for c ollisions be t ween H and H2 have been updated from 
those used in | Shaw et 

3 ( l2005l) . We had used the fi ts given by 
iLe Bourlot et all l ll999l) in Paper 1. lAllers et all l l2005h noted that 
observations of H2 emission from the Or ion Bar suggested tha t 
these rates were two small by nearly 2 dex. IWrathmall et all ( l2007h 
confirmed this and presented a new set of collisional rate coeffi- 
cients using an improved interaction potential and scattering theory. 
These rate coefficients for excited vibrational levels are systemati- 
cally larger, often by 2 to 3 dex. These rates are employed in this 
work. 

The critical density of a level is the density were 
radiative and collisional deexcitation are equally probable 
dOsterbrock & Ferlan j|200 ^. The total radiative rate out of level 
i, Aj, is the sum of the transition probabilities out of the level, 
A/ = £A;/. If the total deexcitation rate [s^'] due to collisions 

Ki 

with species S with density n{S) is n{S)qi = Y, '^{S)qih where q, / 

Ki 

is the collisional deexcitation rate coefficient (units cm-' s^'), then 
the critical density of a molecule in level / colliding with species 5 
is given by 

n(5,/)=A,/g, = £A,, / Y^qij (3) 

Ki I Ki 

The large increase in the H*' - H2 collision rate coefficients for 
vibrationally-excited levels in the new data lowers the critical den- 
sities of most levels by a large amount. Table [2] gives critical den- 
sities for the upper levels of the H2 lines analyzed in this pa- 
per. The rates are evaluated at 1000 K, a typical H2 temperature. 
Each collider considered in our calcu lations is listed. This ta ble 
should be compared with Table 1 of iLe Bourlot et 

3 (I1999') or 

ISternberg & Dalgamo (198^. Our critical densities are substan- 
tially smaller. This means that a Boltzmann level population dis- 
tribution will be established at considerably lower densities than 
was previously thought necessary. 

There is likely to be a substantial uncertainty in all molecular 
collision rates. This is partially because molecular collisions are 
complicated many-body problems. The current H" - H2 rates, while 
sometimes substantially larger tha n the previously publis hed set, do 
not yet include reactive channels JWrathmall et afcoOTi) . The rates 
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Table 2. H2 lines, levels, and critical densities at 1000 K 



WL 


ID 




J 

J up 




h" 


He 




H'' 


H+ 


28.21 


0-0 S(0) 





2 


509.8 


I.UE-hOI 


I.6IE-1-OO 


3.76E-I-00 


5.16E-I-00 


2.36E-02 


17.03 


0-0 S(l) 





3 


1015 


1 54F+02 


2 VOP+O! 






4.14E-01 


12.28 


0-0 S(2) 





4 


1681 


7 492+02 




^ 17F+n2 


A 29F4-02 


1 52F-^^0 


9.66 


0-0 Sf31 





5 


2504 


2 40F+03 


1 43F4-03 


9 Q1F+03 


■X snF4-03 


5 QfiF-^^0 


8.02 


0-0 S(4) 







3474 




fi 47F4-03 


1 28F4-04 


1 SQF+04 


1 27F-I-01 


6.907 


0-0 S(5) 





7 


4586 


1.51E-H04 


2.50E-I-04 


4.45E-I-04 


5.67E+04 


3.02E-H01 


2.223 


1-0 S(0) 


1 


2 


6471 


4.78E+04 


3.61E-I-04 


9.14E-I-04 


1.07E+05 


2.76E-H02 


2.121 


1-0 S(l) 


1 


3 


6951 


3.60E+04 


3.53E-I-04 


8.52E-I-04 


1.07E+05 


2.57E-H02 


2.423 


1-0 Q(3) 


1 


3 


6951 


3.60E-I-04 


3.53E-I-04 


8.52E-I-04 


1.07E-I-05 


2.57E-H02 


2.033 


1-0 S(2) 


1 


4 


7584 


3.16E+04 


5.53E-I-04 


1.15E-I-05 


1.39E-I-05 


2.15E-H02 


1.957 


1-0 S(3) 


1 


5 


8365 


3.05E-H04 


6.77E-I-04 


1.78E-I-05 


2.10E-I-05 


2.17E-H02 


1.891 


1-0 S(4) 


1 


6 


9286 


3.14E-H04 


1.25E-I-05 


2.69E-I-05 


3.21E-I-05 


1.96E-H02 


2.248 


2-1 S(l) 


2 


3 


12550 


1.84E-H04 


5.61E-I-04 


1.25E-I-05 


1.67E-I-05 


3.19E-H02 


1.748 


1-0 S(7) 


1 


9 


12817 


3.81E-H04 


4.44E-I-05 


8.37E-I-05 


1.04E-I-06 


1.35E-H02 



will become larger still, at temperatures where these reactions can 
occur, when this process is included. 

These critical densities can give an indication of the density of 
the molecular gas jjaffe et all200 11. The distribution of level popu- 
lations within H2 ro-vibrationa l levels can be de termined from rel- 
ative emission-line intensities. Ijaffe et all JioOlh find that low v,J 
populations can be fitted as a single Boltzmann distribution cor- 
responding to excitation temperatures --^ 2000 K. In this paper we 
fit the observed H2 line intensities and do not present H2 level- 
excitation diagrams. The results are fully equivalent but have the 
simplification that we are working with the directly-observed quan- 
tity. 

A thermal population distribution will result if the gas density 
is gre ater than the critical densities of the levels involved. I Jaffe et al 
1 I2OOI1) argue that the density in the filaments must be ^ 10* cm ^ 
if the gas is predo minantly molecular. The y based this on a modi- 
fied scaling of the lMandv & Martini ( Il993h semi-c lassical H2 - H2 
collisi on rates. The rate coefficients presented bv lLe Bourlot et all 
( Il999h would suggest a significantly higher density (^10^ cm 
Paper 1 showed that there are other ways to establish a quasi- 
thermal distribution at fairly low densities. Table |2] shows that the 
well-observed 1-0 S(l) line at 2. 121 li-m has a critical density rang- 
ing from 3 X 10^ cm^^, for collisions with H+, to 1.9 x 10^ cm^-', 
for collisions with H2. In practice the critical density will depend on 
whether ions are mixed with the molecular gas. We show below that 
the extra heating and cosmic ray cases have very different mixtures 
of molecules, atoms, and ions, and this may provide a discriminant 
between the energy sources. Note that only collisions with H'* and 

are capable of in ducing a reactive ortho-para transition in H2 
( lDalgarnoetalll973h . 



the equations of statistical equilibrium. The population of the H" 
2p and Is terms determine the excitation temperature 



Tex{2p, Is) 



In 



n{ls)/g{ls) 



n{2p)/g{2p) 



(4) 



where x and k are the excitation potential of 2p and the Boltzmann 
constant respectively. This establishes the sour ce function within 
the L a transition, which we treat following lElitzur & Ferland 
( Il986h . It follows that within the optically thick part of the line 
Sy = Byijex) where Sy is the line source function and By {Tex) is 
the Planck function evaluated at the level excitation temperature. 

Numerical solutions of La transfer in the limits of very large 
optical depths show that the line source function is constant over 
a line width that depends on the total opti cal depth f Adams" 197j; 
iHarringtonl [19731 : iFerland & Netzedll979l) . The results of .Adami 
( I1972I) are adopted. The resulting expression for the full line width 
in the case of large optical depths, t ^ 1, is 



Av = ^VDop 2.3 {ax) 



1/3 



(5) 



where a is the line damping constant and Avoop is the Doppler 
width. With these definitions the photon occupation number within 
the La line is 



riv=Sy 



'(2/ivVc2) = [exp(/iv/;:r„) - 1]" 



(6) 



The rate of induced radiative excitation of H2 by La photons is 
given by 

r> I A ^l-l' Sit A SU —3 —1 ^r^^ 

n,u = niBi Ji,,, = ;j/A„j ^ ^ — = niAuj-qy — cm s . (7) 



'2hv^/c^ g, 



gl 



This is included in the H2 balance equations as a general level ex- 
citation process tShaw et alii2005i) . 



2.5 La pumping of H2 

La fluorescence can be a significant H2 excitation source 
when ions and mo l ecules are mixed as discussed by 
iBlack & van Dishoeckl ( ll987l) and ILupu. France & McCandlissI 
(l2006h . In effect H2 offers a second channel, competing with 
dust, to remove La photons created by collisional excitation 
and recombination. The La energy is converted into internal H2 
excitation energy and eventually dissociation or H2 line emission. 
We determine the full set of H" level populations by solving 



3 THE DENSITY - HEATING PLANE 
3.1 General considerations 

The magnetic field and cosmic-ray density within the filaments are 
not known. The gas density is low enough for the optical [S II] lines 
to be near their low-density limit (Hatchet al 2006), n < 10^ cm^-', 
but the density of the molecular gas is only constrained from 
the form of the H2 level-population distribution jjaffe et aill200lh . 
Given these uncertainties we explore a broad range of non-radiative 
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heating rates and gas densities, compute the thermal and chemical 
state of the gas, and predict the resulting emission-line spectrum to 
compare with observations. Results are presented below as a series 
of contour diagrams showing various predictions. 

The non-radiative energy sources are the dominant heating 
agents for most conditions shown here. For the very lowest val- 
ues of the heating rates the background metagalactic radiation field, 
described above, becomes significant. The resulting emission will 
be unobservably faint in this limit, however, as shown in diagrams 
below and in Paper 1. 

Figure [T] shows the gas kinetic temperature, Figure [2] shows 
the hydrogen molecular fraction, and Figure [3] shows the fraction 
of H in H+ . In all pairs of figures the upper and lower panels show 
the cosmic ray and extra-heating cases respectively. 

An 'effective' ionization parameter, given by the ratio of the 
non-radiative heating to the gas density Unr = rnr/n{ii), character- 
izes physical conditions in these environments. The non-radiative 
heating acts to heat and ionize the gas, while cooling and recombi- 
nation processes are often two-body processes which increase with 
density. The ionization and temperature tend to increase with «„r 
along a diagonal running from low heating and high density, the 
lower right corner, to high heating and low density, the upper left 
corner. The temperature and molecular fractions tend to change 
along diagonals from the lower left to upper right corners. These 
are along lines of roughly constant «,„•. 

The lower right and upper left corners of the plane have the 
most extreme conditions. Gas is cold and molecular in the lower 
left high density, low heating, corner. The kinetic temperature has 
fallen to nearly that of the CMB in this region. A small degree of 
ionization is maintained by background light or cosmic rays. This is 
to be contrasted with gas in the upper left high-heating, low-density, 
corner. That gas is hot, with T ^ 10'* K, and highly ionized. Gas 
in either region would emit little visible/IR light and so would be 
difficult to detect. 



3.2 The H2 emission spectrum 

Paper 1 described the H2 emission spectrum in some detail. Al- 
though the detailed predictions have changed due to the improved 
coUisional rate coefficients used in this paper, the overall trends of 
H2 emission shown in that paper have remained the same. 

The full H2 spectrum is presented below. Here we concentrate 
on the H2 0-0 S(2) A 12.28 (xm line which we will combine with 
optical H I lines to set constraints on our simulations. The log of 
the emissivity Anj (erg cm^^^ s^') of the H2 line is shown in Fig- 
ure |4l H2 emission occurs within regions where the H2 molecular 
fraction, shown in Figure|2l is significant. Other factors besides the 
H2 density determine the line's brightness. 

The emissivity of an H2 line that is excited by thermal colli- 
sions is proportional to 



Log kinetic temperature 



(8) 



where is the density of all colliders and x is the excitation po- 
tential of the upper level of the transition. Higher temperatures and 
H2 densities tend to make the line more emissive. Figures [T] and 
|2]show that «(H2) and T are anti-correlated. From this result, to- 
gether with the form of equation [8] it follows that there will be a 
band of peak emission that runs from the upper right to lower left 
in Figure|4] This is parallel to contours of constant temperature and 
molecular fraction shown in Figures [T] and |2] 

The 11 2. 28 \im line is a pure rotational transition with an level 
excitation energy of 1681 K. The line is excited by two processes 




Log hydrogen density cm' 

Figure 1. The log of the gas kinetic temperature is shown as a function of 
the hydrogen density, the independent variable, and (top panel) the cosmic- 
ray density relative to the Galactic background and (bottom panel) the lead- 
ing term in the extra-heating rate. The temperature is near the CMB in the 
lower right comer and rises along a diagonal from lower right to upper left 
as the heating increases and the gas density decreases. There are discontin- 
uous jumps in the temperature as the gas changes thermal phase. 



in our simulations. The line is emitted efficiently when the gas be- 
comes warm enough for thermal particles to coUisionally excite 
the upper level but not so warm as to dissociate the H2. This pro- 
duces the region of peak emission that tracks the contour giving 
T ~ 2000 K in FigureQ 

The H2 line remains emissive in much cooler regions in the 
cosmic-ray case. Here suprathermal electrons excite the line by a 
two-step mechanism that is analogous to the Solomon process in 
starlight-excited PDRs. The first step is a suprathermal excitation 
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Figure 2. The log of the H2 fraction is shown as a function of the hydrogen 
density and (top panel) the cosmic-ray and (bottom panel) the extra-heating 
rates. Gas is molecular in the lower right corner where the density is high 
and the non-radiative heating rates are low. The gas is nearly fully molecular 
in the lower right corner and becomes increasingly ionized along a diagonal 
running to the upper left. There ai'e discontinuous jumps in the chemistiy as 
the gas changes phases. 



from the ground electronic state X to one of the excited-electronic 
states. A minority of these return to rotation-vibration levels within 
X. Many of these eventually decay to the upper level of the 11 2. 28 
\im line. Other electronic excitations lead to dissociation which 
are then followed by formation by the grain-surface or routes. 
These formation processes populate excited levels within X which 
then lead to emission. Details o f our implementat ion of these H2 
excita tion processes are given in lShaw et all ( |2005|) and lShaw et all 
1 I2OO8I) . 
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Log hydrogen density cm 



Figure 3. The log of the fraction is shown as a function of the hydrogen 
density and (top) the cosmic-ray and (bottom) the extra-heating cases. Gas 
is ionized in the upper left corner of the diagram, where the density is low 
and the non-radiative heating rates are high, and molecular in the lower right 
comer 



3.3 The H I emission spectrum 

We consider the emissivity of H/3 A486lA in detail since we will 
use the Hi 14861A / H2 112.28 |j.m line ratio to constrain our 
simulations. We also consider two H I line ratios that can be used 
to indicate whether the spectrum has been affected by interstellar 
reddening. This is important because the dust content of the fila- 
ments is currently unknown. 

Figure [5] shows the log of the emissivity Ak j of H/3 A4861 A, 
one of the brightest lines in the optical spectrum. In a photoionized 
environment H I lines form by radiative recombination of + H+ . 
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Log 4jty(H2 12.28 nm) (erg cm"^ s"^) 




1 2 3 4 5 6 



Log hydrogen density cm' 

Figure 4. The log of the emissivity Anj (erg cm"^ s"') of the Hi 0-0 S(2) 
12.28 jjm hne is shown as a function of cloud parameters. This figure is 
to be contrasted with the next figure showing the H I hne emissivity over a 
similar range of parameters. 

H I lines also form by recombination if the gas is collisionally 
heated to a temperature so warm that nearly all H is ionized. The 
emissivity of an H I recombination line is proportional to 

4;ry(Hi)oc„(H+))ier"°'^ (9) 

where T is the gas kinetic temperature and the inverse temperature 
dependence is a fit to recombination coefficients in the neighbor- 
hood of 10^ K. H I recombination lines are emitted most efficiently 
by gas that is dense, ionized, and cool. 

H I lines can also form by collisional excitation in regions 
where H is mostly atomic. The gas must be warm enough to excite 
the upper levels of the optical or IR H I lines, requiring a tempera- 



Log 47iy(Hp) (erg cm'^ s'^) 




1 2 3 4 5 6 



Log hydrogen density cm" 

Figure 5. The log of the emissivity Anj (erg cm"^ s"') of H/5 /l486lA is 
shown as a function of cloud parameters. This figure is to be contrasted with 
those in Paper 1 which show the H2 line emissivity over a similar range of 
parameters. 

ture T ^ 4000 K, although collisional excitation is important at low 
temperatures in the cosmic ray case when suprathermal secondary 
electrons are present. The emissivity of a collisionally excited H I 
line is proportional to «(H*')%/(r) where n{Yf') is the atomic hy- 
drogen density, mostly in the ground state, and % is the density of 
the colliding species. The function f{T) is the Boltzmann factor of 
the upper level of the H I line for the case of excitation by thermal 
particles and is a constant for excitation by secondary electrons. 
Collisionally excited H I lines are emitted most efficiently by gas 
that is dense, atomic, and warm. 

The extra heat case is relatively simple since only thermal pro- 
cesses affect the excitation and ionization of the gas. The lower 
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panel of Figure |5] shows that the hne has a ridge of relatively high 
emissivity which corresponds to moderate ionization and tempera- 
ture. The H I lines along this ridge are mainly produced by colli- 
sional excitation with a contribution from recombination following 
coUisional ionization. The emission falls off dramatically when the 
gas becomes molecular in the lower right corner or very hot in the 
upper left comer. 

Figure|5]shows the H/3 emissivity in the cosmic-ray case. Col- 
lisional excitation is also a dominant contributor to the line. The 
ridge of peak emission occurs for the same reason as in the ex- 
tra heating case. Significant H I emission occurs in the lower right 
comer since the dissociative cosmic rays prevent the gas from be- 
coming fully molecular even when it is quite cold. The population 
of secondary electrons produces a significant coUisional excitation 
contribution to the H I lines across most of the lower right half of 
Figure [5] 

In both heating cases the emissivity tends to rise along a diag- 
onal from the lower left to upper right. This corresponds to rising 
iieiip at nearly constant T. This figure clearly indicates that power- 
ful selection effects operate in this environment. As we show below, 
the emissivity of the H I lines is substantially higher than produced 
in the pure-recombination case due to the coUisional contribution. 
This affects the use of H I lines as indicators or the gas of ionized 
gas. 

The H I lines are often used to measure the amount of in- 
terstellar extinction. It is unusual for the H I recombination spec- 
trum to deviate very far from Case B for moderate densities and 
optic al depths in the Lyman continuum in a photoionized cloud 
jOst erbrock & Ferland 2006). This is because the relative H I line 
intensities are mainly determined by the transition probabilities. 
These determine which lines are emitted as the electrons cascade 
to ground after capture from the continuum. The result is that the 
spectrum depends mainly on atomic constants and less so on the 
physical conditions in the gas. 

The ratio of the intensity of Ha A 6563 A relative to 
H/3 14861 A, two of the strongest lines in the optical spectrum and 
a possible reddening indicator, is shown in Figure |6l The Ha/Hj3 
intensity ratio is ~ 2.8 for Ca se B in low-density photoionized gas 
dOsterbrock & Ferland |200^ . The predicted ratio is significantly 
larger than this in most regions where the H I lines have a large 
emissivity. 

It is possible to measure intensities of H I lines that originate 
in a common upper atomic level with spectra that cover both the 
optical and infrared. The Pa 1 1 .87|i.m and H/3 A4861 A lines, with 
a common n = 4 upper configuration, is an example. Each of these 
lines is actually a multiplet with 4 nl terms within the upper level. 
At conventional spectroscopic resolution the lines within the mul- 
tiplet appears as a single H I line. The transition probability of this 
multiplet depends on the distribution of nl populations with the up- 
per terms, which in turn depend on the gas density ( Pengellv 1964i). 
Such line ratios are not expected to depend on the physical condi- 
tions as much as lines that originate in different configurations. 

Figure |7] shows the computed Pa H.87iJ.m / H/3 A486lA in- 
tensity ratio. It is not constant because of the changing populations 
of nl terms within the n = 4 configuration of H". At high densities 
the / terms will be populated according to their statistical weight 
and the level populations are said to be well /-mixed. The line in- 
tensity ratio is then a constant that is determined by the ratio of 
transition probabilities and photon energies. Cloudy treated n> 2 
configurations in this well /-mixed limit in versions before C08, the 
version described in this paper. As described above, these calcula- 
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Figure 6. The ratio of emissivities of H I Ha A6563 to H/3 A4861 is shown 
as a function of cloud parameters. The ratio is ~ 2.8 under the Case B 
conditions expected for photoionized nebulae. It is significantly larger for 
most parameters shown here because of coUisional contributions to the line. 



tions use an improved H model which fully resolves the nl terms 
within a configuration. 

When the density is too low to collisionally mix the nl terms 
their populations tend to peak at smaller / since recombinations 
from the continuum and suprathermal excitation from ground tend 
to populate low-/ levels. It is only when the density becomes large 
enough for /-changing collisions to mix the / levels that the well /- 
mixed fimit is reached. This makes the Pa A 1. 87 nm/ H/3 14861 A 
intensity ratio depend on density. 

Figure |7] shows that the Pa / H/3 intensity ratio does not vary 
over a broad range. The extra-heating case, shown in the lower 
panel, has a Pa / H/3 ratio in the range of 0.2 — 0.3, below the Case 
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Log hydrogen density cm' 

Figure 7. The H I Pa A1.87|am to Hj3 A4861 emissivity ratio is sliown as 
a function of cloud parameters. These lines have a common upper n = A 
configuration so changes in the relative intensities are due to changes in the 
nl populations within the n = 4 configuration. The Case B value is 0.34. 



B ratio of 0.34 at 10'* K and low densities. The cosmic ray case has 
a Poc / Hj3 intensity ratio that varies between ~ 0.2 — 0.4, within 
a factor of two of the Case B ratio. The range is larger because 
suprathermal electrons excite the atom for nearly all temperatures. 
These non- thermal excitations are mainly \.o np terms and create a 
distribution of nl populations that differs from the recombination 
case. This line ratio remains a good reddening indicator because of 
the relatively modest range in the predicted values when combined 
with the very wide wavelength separation of the two lines. 



Table 3. Physical conditions for the two cases at a range of kinetic temper- 
atures. 



Species 


r = 10^ K 


T ~ 700 K 


r = 10'' K 


n(,(Heat, cm"^) 


0.013 


0.0089 


0.59 


>J(,(CR, cm^^) 


0.92 


9.3 


155 


n(H2)(Heat, cm^^) 


952 


940 




n(H2)(CR, cm-3) 


232 


16 




n(H'')(Heat, cm^^) 


48 


60 


999 


n(H'')(CR, cm-3) 


768 


975 


863 


n(H+)(Heat, cm^^) 






0.0004 


«(H+)(CR, cm-3) 


0.80 


8.4 


137 



3.4 Chemical and ionization state of the gas 

The two non-radiative heating cases, while having similar overall 
trends, have very different detailed properties, which we examine 
next. Figure [8] shows predictions for a vertical line in Figures [T]- 
[3]corresponding to ne = 10^ cm^^ and varying non-radiative heat- 
ing rates. The independent axis gives the extra heating rate pre- 
coefficient and the cosmic-ray density relative to the Galactic back- 
ground. These parameters were chosen to cross the region where 
the physical conditions in Figures [T]|2] and [3] change abruptly. The 
left panels show the extra heating case and the right panels the cos- 
mic ray case. The second from the bottom panels of Figure[8]show 
the kinetic temperature. The discontinuous jump in temperature, 
the regions where the contouring in Figures[T]-[3]become ragged, 
is clear. The lowest panels show the molecular, atomic, and ion- 
ized fractions of hydrogen. The upper right panel of Figure[8]shows 
some heating and excitation efficiencies for the cosmic-ray case. In 
all cases the temperature and ionization of the gas increases to the 
right as non-radiative energy sources become more important. 

Both non-radiative heating cases have discontinuous jumps in 
temperature. These are due to inflection points in the gas cooling 
function that occur when more than one solution is possible. These 
stability issues are discussed further below. We first focus on the 
changes in the ionization of the gas shown in the lower two panels. 

Table [3] compares physical conditions at three temperatures. 
The low and high temperatures of 100 K and lO'* K are near the 
extremes of the regions which produce the molecular and low- 
ionization emission. The mid-temperature was evaluated as close 
to lO"' K as possible. This is a typical temperature for H2- emitting 
gas, as shown in Paper 1. No stable thermal phase exists at exactly 
this temperature so the conditions at 700 K, the warmest stable re- 
gion below 10^ K, IS given. 

In the extra-heating case thermal collisions are the only ion- 
ization source. The result is that at low temperatures the gas is pre- 
dominantly molecular with a trace of H*'. As the heating rate and 
temperature increase there is an abrupt change in the chemical state 
when the kinetic temperature approaches the dissociation potential 
of H2. The dissociation of H2 causes the particle density to increase 
and the mean molecular weight to decrease. Both cause the colli- 
sional rates, excitations, cooling, and ionization, to increase. The 
result is a positive feed-back process that causes an abrupt phase 
transition from H2 to H*' . This is the reason that contours overlap 
in the lower panel of Figures [T]-[3] 

This behavior is in sharp contrast with the cosmic-ray case. 
The effects of relativistic particles on a predominantly thermal gas 
has been well documented in a number of papers starting with 
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Figure 8. This shows the thermal and physical state along a vertical line at nw = 10 cm^ in the previous contour plots. The left panels are the extra-heating 
case while the right panels are the cosmic-ray case. The non-radiative rates are the independent axis in each panel. The range in both was adjusted so that the 
phase transition, where many physical quantities change abruptly, occurs near the middle of the independent axis. The ionization and temperature increase to 
the right as the non-radiative heating rates increase. The top right panel shows the cosmic ray heating, ionization, and line-excitation efficiencies. Cosmic rays 
ionize and excite predominantly neutral gas and heat ionized gas. The abrupt change in temperature, shown in the second to bottom pair of panels, is more 
extreme in the extra-heating case. The two lower panels show the physical state of hydrogen. There is a mix of atomic and molecular gas in the cosmic-ray 
case due to the ionization and dissociation that they produce in cold neutral gas. 



ISpitzer & Tomaskol 1 1968h and most recently by IXu & McCravl 
(Il99lh and lDalgamo et ^19991) . As shown in the upper right panel 
of Figure[8l cosmic rays excite, heat, and ionize the gas. The curves 
marked 'Ionization' and 'excitation' show the ionization and ex- 
citation rate (s^') but have been divided by the cosmic-ray-to- 
background ratio to remove the effects of increasing cosmic-ray 
densities. The curve marked 'heating' has been similarly scaled and 



further multiplied by 10^ to facilitate plotting. The fraction of the 
cosmic-ray energy that goes into each process depends mainly on 
the electron fraction, yielny^, the ratio of the thermal electron to total 
hydrogen densities. This electron fraction increases as the cosmic- 
ray rate increases. 

For low electron fractions a cosmic ray produces a first gen- 
eration secondary electron that is more likely to strike atoms or 
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molecules causing secondary excitations or ionizations. Radiation 
produced by line excitation or recombination following ionization 
will be absorbed by dust which then reradiates the energy in the 
FIR. Relatively little of the cosmic-ray energy goes into directly 
heating the gas. As the electron fraction increases, moving to the 
right in Figure [8] more of the cosmic-ray energy goes into heating 
rather than exciting or ionizing the gas. This is because for larger 
electron fractions the secondary electrons have a greater probabil- 
ity of undergoing an elastic collision with a thermal electron. This 
adds to the thermal energy of the free electrons and so heats the 
gas. 

The effect is to produce a more gentle change in the ioniza- 
tion of the gas as the cosmic-ray rate is increased. Ionization and 
dissociation are mainly caused by non-thermal secondary particles. 
Significant levels of dissociation or ionization exist at temperatures 
that are too low to produce these effects by thermal collisions. The 
lower two panels of Figure [8] show how the H2 and H" fractions 
change . In the cosmic-ray case they change continuously until 
the thermal instability point is reached. Although a jump still oc- 
curs, the change in temperature is mitigated by the more continuous 
change in ionization and molecular fractions. This is to be con- 
trasted with the extra-heating case where the gas is almost entirely 
molecular or atomic with little mixing of the two phases. 

Some specific values are given in Table[3l At the lowest tem- 
perature, 100 K, nearly all H is molecular in the extra-heating case, 
while in the cosmic-ray case substantial amounts of H" are present. 
There are even significant amounts of H+ at this temperature be- 
cause of the ionization produced by the cosmic ray secondaries. At 
the highest temperature shown, 10* K, nearly all H is atomic in the 
extra heating case while a substantial amount of H+ is present in 
the cosmic-ray case. 

A substantial population of H'* and H+ is mixed with H2 in 
the cosmic-ray heated gas. Molecular regions become far warmer 
and have far fewer H+ ions mixed with H2 in the extra-heating 
case. This leads to the most important distinctions between the two 
cases. Ions have larger collision cross sections so are more active 
in altering the level populations of H2, as shown in Table |2] The 
H2 populations will be different as a result. Further, H'^ and H+ 
can underdo orth-para exchange collisions with H2. This process 
will be far more rapid in the cosmic-ray case and will be another 
distinction between the two cases. 



3.5 Thermal state of the gas 

Figure [8] shows that both cases have temperatures that change dis- 
continuously. The origin of these jumps is described here. We con- 
centrate on the extra-heating case where the effects are the largest. 

Figure |9] shows gas coolants as the extra heating rate is var- 
ied across Figure [8] The coolants leftward of the discontinuous 
break are classical atomic and molecular 'PDR' coolants. A Galac- 



tic PDR is the H" re 



'ion adjacent to galactic regions of star for- 



mation dTielensllIool . These lines can be detected by current and 
planned IR instrumentation. 

The gas abruptly changes from H*^ to H+ at the discontinu- 
ity where the temperature jumps from 'warm' 10"* K) to 'hot' 
(~ 10* K). The coolants on the hot side are mainly emission lines 
of atoms and ions of the more common elements. The strongest 
coolant, marked 'H I lines', is the set of Lyman lines that are col- 
lisionally excited from the ground state. Although their intrinsic 
intensity is large they do not escape the cloud because of the large 
H I line optical depths and the presence of dust. These lines are 
absorbed by dust as the photons undergo multiple scattering. This 




Figure 9. This identifies the most important cooling transitions for the 
extra-heating case shown in Figure [8] The cooling shifts from the IR when 
the gas is cold and molecular, regions with log//o < —23, into the optical 
and UV when the gas is warm and ionized. 



is an example of a process that cools the gas by initially convert- 
ing free electron kinetic energy into line emission which is then 
absorbed by, and heats, the dust. The energy eventually escapes as 
reprocessed FIR emission. Note that the gas and dust temperatures 
are not well coupled for the low densities found in these filaments. 
The dust is generally considerably cooler than the gas. 

The behavior shown in the two previous figures, with discon- 
tinuous jumps in the physical conditions, is due to well-known ther- 
mal instabilities in interstellar clouds. Figure [Tol shows the calcu- 
lated cooling rate for a unit density and, again, the extra heating 
case. The heating rate is varied and the temperature determined by 
balancing heating and cooling. We plot the cooling as a function of 
temperature rather than the heating rate to make it easier to com- 
pare these results with previous calculations. 

The overall shape of t he cooling function i s discu ssed, for 
instance, in the review by iDalgarno & McCravl ( Il972h . At low 
temperatures cooling is mainly by molecules and low-lying lev- 
els within ground terms of atoms and first ions of the abundant 
elements. At low temperatures cooling usually involves changes in 
levels which have excitation energies of ^ 10"' K. The cooling in- 
creases as the Boltzmann factor for excitation increases and reaches 
a peak at ~ 10^ K where the Boltzmann factor reaches unity. The 
cooling then declines as T increases further due to the r ^ ' depen- 
dence of the Coulomb focusing factor. Starting at roughly 5000 K 
lines involving changes in term become energetically accessible, 
increasing the cooling, and producing a peak at --^ 10* K. The de- 
cline after the peak is again due to the r^'/^ dependence of the 
collision rate coefficient when the Boltzmaim factor has reached 
unity. The similar peak at ~ 10^ K is due to transitions involving 
changes in electronic configuration in ions of second and third-row 
elements. 

The form of the cooling curve is affected by the presence of 
dust. Many of the most important coolants in a dust-free mixture, 
such as iron, calcium, silicon, and others, are strongly depleted 
when dust forms. The result is that cooling is less efficient due to 
the loss of these gas-phase coolants. This is the major reason that 
the cooling curve differs from the solar-abundance case. 

The condition for thermal stability of a constant-density gas is 
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Figure 10. The cooling function for gas with unit density and a range of 
extra-heating rates. The heating is varied and the equilibrium temperature 
determined. The derived gas kinetic temperature is used as the independent 
variable to better compare with previous calculations and so that thermal 
stability can be judged. The arrows indicate regions where the gas will un- 
dergo a discontinuous jump in temperature to avoid thermally instable re- 
gions. A gas that was originally cold and molecular would follow the curve 
moving from left to the right while initially hot gas would move from right 
to left. 
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( lFieldlll965h . Here C and H are the cooling and heating rates, erg 
cm^^ s^^, and p indicates that the derivative is at constant density. 
Small perturbations in the temperature cause the temperature to 
continue to change in the direction of the perturbation if the deriva- 
tive is negative. Our hypothesized non-radiative heating processes 
have no explicit temperature dependence so the partial derivative 
with respect to temperature simplifies to the derivative of the cool- 
ing. The portions of the cooling curve shown in Figure [TO] with 
negative slope are thermally unstable. 

Thermal instabilities cause the gas to have a memory of its 
history. The present state of the gas will depend on whether it 
approaches thermal equilibrium from the high or low-temperature 
state. One possibility is that the filaments c ooled down from the sur- 
rounding hot-ionized plasma tRevaz et a I2OO&) . They would have 
reached their current state after moving from right to left in Fig- 
ure [To] In this case the gas would follow the leftward-pointing ar- 
rows when it passed the unstable regions. Filaments originating as 
cold molecular gas, perhaps in the interstellar medium of the cen- 
tral galaxy, would move from left to right and would follow the 
rightward-pointing arrows. Different regions of the cooling curve 
are reached by a heating or cooling gas. This difference could pro- 
vide a signature of the history of the gas. 



3.6 The grain/molecule inventory and the history of the gas 

The dust content and molecular inventory provides a constraint on 
the history of the filaments. If the filaments formed from the sur- 
rounding hot X-ray plasma they would likely be dust free. Proper- 
ties of du st-free clouds within a galaxy cluster were examined by 
iFerland et al ( 1994) and Ferland et al (2002). If the gas originated 
in the interstellar medium of the central galaxy it would most likely 
contain dust, as we have assumed so far. 

The timescales required to form the observed molecules pro- 



vide a clue to the origin of the gas. The slowest step is the for- 
mation of H2. There are no direct formation processes because the 
homonuclear molecule has no permanent dipole moment. In dusty 
environments H2 forms by cata l ysis o n grain surfaces. We have 
adopted the rate derived by ljural jl975l) for the galactic ISM. The 
H2-formation timescale is then 



10° 



lO^cm-3 Al 



years 



(11) 



where «(H'^) is the atomic hydrogen density and A^i,,, / A/sm is the 
filament dust to gas ratio relative to the ISM value. 

If grains are not pres ent then H2 will m ainly form by asso- 
ciative detachment of 1 Ferland et 3|l99i). This process is lim- 
ited by the rate of the slowest step, radiative association to form 
H^. This has a rate coefficient given by r = 8 x lO^^^t'^'^^ n,, s^' 
in the neighborhood of 10"^ K. Carbon is the most abundant electron 
donor in H'^ regions where this process is fast. The timescale can 
be posed in terms of the hydrogen density by assuming a solar C/H 
ratio and that all C is singly ionized. The H2 formation timescale in 
a dust-free environment is then 
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where the density, temperature, and C/H abundance are given in 
units of 10^ cm^^, IQ-' K, and solar, respectively. 

The dust-free H2 formation timescales are ~ 2.5 dex longer 
than those in a dusty environme nt. Unfortu nately the ages of the 
filaments are poorly constrained. iHatch et aU2006.) find that a fila- 
ment in the Perseus cluster with a projected length of ~ 25 kpc and 
a velocity spread of ~ 400 km s^' . The corresponding expansion 
age is > 7 X 10^ yr. This is comfortably longer than the formation 
time for a dusty gas. H2 could only form in the dust-free case if the 
density is considerably higher than expected from typical pressures. 
This is possible but unlikely. We return to this question below. 

The stable thermal solutions presented in the remainder of this 
paper are limited to those on the rising-temperature branch of the 
cooling curve. This, together with our assumed dust content, im- 
plicitly assumes that the gas was originally cold and dusty. In prac- 
tice this means that the temperature solver determines the initial 
temperature with a search that starts at low temperatures and in- 
creases T to reach thermal equilibrium. 



4 PROPERTIES OF CONSTANT-PRESSURE CLOUDS 
4.1 Equation of state 

The equation of state, the relationship between gas density and the 
temperature or pressure, is completely unknown for the filaments. 
Possible contributors to the total pressure include gas, turbulent, 
magnetic, cosmic ray, and radiation pressures. The gross structure 
of the filaments, forming large lines or arcs, is reminiscent of mag- 
netic phenomena like coronal loops. This suggests that magnetic 
fields guide the gas morphology with the field lines lying along the 
arc. Gas would then be free to flow along field lines, the long axis 
of the filaments, but not across them. 

We assume that the gas pressure within the filaments is 
the same as the gas pre ssure in the surrounding hot gas, nT = 
W^-^ cm"3 K, found by I Sanders & FabianI ilOOll) in the Perseus 
cluster. The magnetic field must then be strong enough to guide the 
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cool material and so maintain the linear geometry despite any mo- 
tion with respect to the surrounding hot gas. A magnetic field of 
~ 100/iG has an energy density equivalent to a thermal gas with 
nT = 10^'^ cm^-' K so a field stronger than this is needed. Little is 
known about the filament geometry at the sub-arcsec {8r ^ 35 pc 
at Perseus) level and nothing is known about the magnetic field 
strength within the filaments. Magnetic confinement does provide 
a plausible explanation for the observed geometry however 

Figure [TT] shows the gas pressure in conventional ISM 
unit {P/k = nJcm^-^K) for the two cases. If individual sub- 
components that make up the filaments have constant pressure then 
they will lie along one of these contour lines. 

The dynamic range in Figure[TT]is large. The solid line in Fig- 
urel 1 2lcorresponds to our preferred pressure of riT = lO*"'^ cm^^ K. 
The dashed lines in the figure show pressures 0.5 dex to either side 
of this value. As we show next, regions where the contours exceed 
a 45 degree angle, and where the ±0.5 dex contours nearly overlap, 
are thermally unstable. 



4.2 Thermal stability 

The criterion for thermal stability in a constant-pressure gas is 
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jFieldl 19651) . where the equality assumes constant mean molecular 
mass. The first term on the right is the same as that considered in 
the constant-density case discussed in the Section [33] Both of the 
heating sources considered here have linear density dependencies 
and no explicit dependence on temperature. At constant pressure 
this will carry over into a heating rate that is inversely proportional 
to the temperature so the second term is non-zero. 

The upper panel of Figure[T3]shows the gas cooling rate along 
the isobaric log nT = 6.5 contour in Figure[T2l Both heating cases 
are drawn in each panel. The figure shows the cooling as a function 
of temperature rather than the heating rate so that the gas stabil- 
ity can be more easily judged. The detailed cooling properties of 
the two sources of non-radiative heating are quite different, as can 
be seen from the physical properties shown in Figure [8] and listed 
in Table [3] The line marked 'temperature^'' shows the approxi- 
mate form of the second term in equation [T3] Regions where the 
slope of the cooling function is steeper than this are unstable. The 
lower panel of Figure[T3] shows the product of the cooling and the 
temperature. Expressed this way, regions with negative slope are 
unstable. Two distinct phases, corresponding to cold molecular and 
warm atomic/ionized, exist. 

The upper axis in Figure [T3] and the figures that follow, gives 
the particle density along the isobaric line while the lower axis 
gives the kinetic temperature. Temperature and density are related 
since the product nT is constant. Note that the density n is the to- 
tal particle density not the hydrogen density. For molecular regions 
the total particle density is about half the hydrogen density while 
in ionized regions it is about twice due to the presence of free elec- 
trons. 

Figure [T4l shows the cooling time for the cases shown in Fig- 
ure [T3] The cooling time will determine how quickly the gas can 
respond to changes in the environment. It also determines how 
quickly thermally unstable gas will heat or cool and reach stable 
regions of the cooling curve. These times are short compared with 
timescales over which the galaxy cluster changes and are far shorter 
than the H2 formation timescales mentioned above. 



Log nT 




Log hydrogen density cm 



Figure 11. The log of the gas pressure, expressed as P/k = nT(Kcm^^), 
is shown as a function of the hydrogen density and (top panel) the cosmic- 
ray rate relative to the Galactic background and (bottom panel) the extra- 
heating rate. For compaiison inner regions of the Perseus cluster have nT « 
10*"^ cm-3 K. 



4.3 Spectrum emitted by a homogeneous cloud 

The full emitted spectrum is computed for each point along the 
isobaric line in Figure[T2] Emissivities 471 j of a few representative 
lines are shown in Figure [Ts] The luminosity of a line will be the 
integral of 4;r j over the volume of the emitting region. 

Most of the lines shown in Figure[T5]will be optically thin for 
any reasonable column density. In this case their emissivity has no 
explicit dependence on column density. If we can neglect internal 
reddening, a good approximation for the IR and radio lines but an 
additional uncertainty for optical or UV transitions, then the line 



© RAS, MNRAS OOO.dUlT] 



Collisional heating of cool cluster filaments 1 5 



Log nT 



Log particle density cm" 
5 4 3 






2 3 
Log temperature K 

Log particle density cm"'' 
5 4 3 



10-' 



o 

I 10" 

o 
o 
O 



10' 



10-' 



Extra Heating 




Log hydrogen density cm" 



2 3 
Log temperature K 

Figure 13. The solid lines in tlie upper panel give the volume cooling rates 
for the two cases as a function of the temperature (the lower independent 
axis) and particle density (he upper axis) along the P/k = nT = 10*^ cm"^ 
K isobaric contour of the previous figure. Regions where the slope of the 
cooling curve is steeper than the dashed line mai'ked temperature^ ' in the 
upper panel are thermally unstable. The lower panel shows the product of 
the cooling and the temperature. Regions where this product has a negative 
slope are unstable. 



Figure 12. The solid contour corresponds to the preferred gas pressure 
P/k = nT = lO*'^ cm^' K. The dotted lines correspond to gas pressures 
0.5 dex above and below this value. Regions where the contours have a 
steeper than 45 degree slope, where the contours nearly overlap, are ther- 
mally unstable. 



luminosity is the integral of the emissivity over the volume con- 
taining gas. This is simple since the luminosity of an optically thin 
line is simply proportional to the total volume of emitting gas rather 
than on details of how the gas is arranged. 

The CO lines are more problematic since lower-y transitions 
are generally optically thick in the Galactic ISM. The emissivity 
of the CO J = 3 ^ 2 XS63[im tra nsition is shown ra ther than 
the lower-y lines normally observed i lSalome et a3l2006h . since it 
is more likely to be optically thin, if the column density of CO 
is large enough to thermalize an optically-thick line then the lu- 



minosity will be set by the gas kinetic temperature and the cloud 
area rather than on the emissivity. It is by no means certain that the 
A 863 vim will be optically thin in the filaments and we shall return 
to the CO line spectrum below. 

Figure [15] shows that different lines form at different tem- 
peratures and densities because, as Figure [8] shows, atoms and 
molecules exist at different regions with little overlap. 

Paper 1 showed that H2 lines in the filaments are much 
stronger relative to H I recombination lines than is found in nor- 
mal star-forming regions like Orion. H I recombination lines will 
have a peak emissivity that is weighted towards the coolest and 
densest regions where H+ is present, as shown in Equation |9] and 
Figure|5] This is indeed the case in the extra-heating case as shown 
in the lower panel of Figure[T5] where we also see that the H I and 
H2 lines form in gas with distinctly different densities and temper- 
atures. 

The extra-heating case, shown in the lower panel of Figurell5l 
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Figure 14. The cooling timescales for the two cases shown in the previous 
figure is shown. Gas in unstable regions of the previous figure will move to 
a stable portion of the cooling curve on this timescale. 

is simplest and we consider it first. Lines excited by thermal colli- 
sions tend to form in warmer gas because of the exponential Boltz- 
mann factor. Low-excitation H2 lines such as H2.28p.m form in 
cooler gas than the H2 A2. 12^i.m, which has an upper level with an 
excitation potential of --^ 7000 K, in the extra-heating case. 

The emissivities of the Ha and H2 lines have two local peaks 
in the cosmic-ray case, as shown in the upper panel of Figure [TS] 
The lower-temperature peak is due to the direct excitation of H I 
and H2 lines by secondary electrons. Excitation to H2 electronic 
levels that then decay into excited states of the H2 gr ound state, 
a process analogous to the photon Solomon process JShaw et^ 
[2005 ). contributes to the H2 emission. The higher-temperature peak 
occurs when cosmic rays heat the gas to a warm enough tempera- 
ture to excite H2 transitions with thermal collisions. Eventually, in 
the highest-r regions of the figure, the gas becomes warm enough 
to be predominantly ionized and the H I lines begin to form by re- 
combination. 

The hydrogen emission-line spectrum is predicted to be differ- 
ent from Case B relative intensities in regions where the lines form 
by coUisional excitation of H" rather than by recombination of H+. 
Figure[76lshows the predicted intensities of the brighter optical and 
IR lines relative to Ha. There are two reasons for differences from 
simple Case B, the first being the large range in kinetic temperature 
along the isobaric line. The standard Case B spectrum is most often 
quoted for T = \(f K, appropriate for a photoionized cloud. The 
fact that coUisional processes excite the lines in this environment 
cause further deviations from Case B. 

Table H compares the H I spectrum for two points along the 
constant-«r line with Case B predictions. Little H I emission oc- 
curs for the extra-heating case at T ~ 300 K so only the cosmic-ray 
case is given. The extra heating and cosmic ray predictions agree 
at r ~ 2 X 10^ K because the gas is warm enough to collisionally 
excite the H I lines. The emissivity 4;r j is --^ 1 dex brighter than 
a pure recombination H I spectrum because of the contribution of 
coUisional excitation of H^ . Note that the quantities given in Table 
|4]are the emissivity A nj and not the emission coeffi cient A7t j/n^tip 
given in, for instance, lOsterbrock & Ferlandl l l200^ The relative in- 
tensities of the H I lines are also different in the two cases. The H I 
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Figure 15. The emissivities, the emission per unit volume, are shown for 
several emission lines along the isobaric line corresponding to P/k = nT = 
10*^ cm"^^ K. The total pailicle density and temperature are indicated as 
the independent axes. The upper panel is the cosmic ray case and the lower 
panel is for extra heating. 



lines at r ~ 300 K are mainly excited by non-thermal secondary 
electrons with collision rates that are proportional to the oscilla- 
tor strength of the Lyman line. The emission at T ~ 2 x 10^ K is 
mainly produced by collisions with thermal electrons whose rates 
are dominated by the near-threshold cross section. Case B applies 
when the lines form by recombination, which is not the case here. 

Table |4] shows that, although the H I lines have a far greater 
emissivity than is given by Case B, the relative intensities are not 
too dissimilar. This is because the relative intensities of the H I 
spectrum are most strongly affected by the transition probabili- 
ties that determine how a highly-excited electron decays to ground 
rather than detail detains of how the electron became excited. 

4.4 Allowed range of cloud density and temperature 

The observed spectrum has a wide range of molecular and atomic 
emission, requiring that cloudlets with a range of densities and tem- 
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Figure 16. The intensities of several strong optical and IR H I lines are 
shown for points along the isobaric line coiresponding to lognr = 6.5. The 
total particle density and temperature are indicated as the independent axes. 



Table 4. Predicted H I emission spectrum 



3.0 



Log particle density cm 
2.5 2.0 



1.5 





CR 


CR,H 


CaB 


Heat 


CR 


logMH 


4 


2.2 


2 


integrated 


integrated 


logr 


2.5 


4.3 


4.3 






4;r;'(Ha) 


-19.85 


-19.96 


-21.01 


-20.029 


-20.171 


H I 3798A 


0.0132 


0.00076 


0.0197 


6.41(-3) 


0.0109 


H I 3835A 


0.0181 


0.0112 


0.0272 


9.38(-3) 


0.0152 


H I 3889A 


0.0261 


0.0172 


0.0391 


0.0145 


0.0224 


H I 3970A 


0.0400 


0.0286 


0.0591 


0.0244 


0.0352 


HI4102A 


0.0665 


0.0531 


0.0962 


0.0459 


0.0608 


H I 4340A 


0.125 


0.116 


0.173 


0.102 


0.120 


Hl486lA 


0.288 


0.254 


0.364 


0.227 


0.271 


H I 6563A 


1.000 


1.000 


1.000 


1.000 


1.000 


Hi 1.875|am 


0.121 


0.0815 


0.104 


0.0839 


0.0992 


Hi 2.166 |jm 


0.0082 


0.0053 


0.085 


4.84(-3) 


6.35(-3) 



peratures occur. Only some of the solutions shown in Figure [Hwill 
exist in a constant-pressure filament. The lowest temperature is set 
by the CMB temperature, the lowest temperature that molecules 
and grains will have. This low-temperature limit sets a high-density 
limit since the product nT = 10^'^ Kcm^-* is constant. Gas denser 
than ncMB ~ 10^ K will remain at the CMB temperature, be over- 
pressurized relative to its environment, and would expand. We do 
not consider gas denser than ncMB • 

Thermal stability (Section 14.21 and Figure [T3] l sets a low- 
density, high-temperature, limit if the system is time steady. Gas on 
unstable parts of the cooling curve will move to regions of thermal 
stability on the cooling timescale of the region. We only consider 
thermally stable solutions here and come back to discuss this point 
below. 

Figure [TTIshows the emissivities of some strong lines that are 
produced near the low-density, high-temperature, end of the cloud 
distribution. We concentrate on the optical spectrum in this Figure 
since these lines form in warm gas and because their relative inten- 
sities can be measured with greater precision owing to the fact that 
this spectral range can be observed with a single entrance slit. The 
grey band in the upper right of each panel indicates regions which 
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Figure 17. The emissivities of several of the prominent optical atomic and 
low-ionization lines seen in the filaments are shown. The lines include lO I] 
A6300, lO II] A3727, IN I] A5199, and [N II] A6584. The shaded rectangle 
indicates thermally unstable regions. 



are thermally unstable. Most of the lines shown in the Figure have 
peak emissivities that occur in stable regions. The predicted inten- 
sities of these lines will not depend greatly on the precise high-r 
cutoff of the cloud distribution. 

Figure [Ts] shows some fainter lines which are currently unob- 
served in the Horseshoe filament in the Perseus cluster which we 
model in detail below. Some, for instance [O III] A5007A, will be 
very sensitive to the thermal stability cutoff since the peak emis- 
sion occurs in unstable regions. We show below that these lines can 
discriminate between the two non-radiative heating processes. 



5 EMISSION FROM A DISTRIBUTION OF CLOUDS 

5.1 The need for clouds with a range of properties 

The spectrum shows that emission from both molecula r and ionized 
gas comes from gas that is nearly spatially coincident dSalome et all 
.2008b) . Physically this requires that gas with a variety of densities 
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Figure 18. Tlie emissivities of several optical atomic and low-ionization 
lines that are either sensitive to the thermal stability cutoff or which might 
discriminate between the two non-radiative heating models we discuss. The 
lines include [O III] A5007, [Ne III] A3869, and He I A5876. Ha is shown 
for reference. The shaded rectangle indicates themially unstable regions. 



contributes to the net emission detected by the entrance aperture of 
a spectrometer. 

No one choice of cloud parameters can reproduce the full 
range of molecular, atom, and ionic emission. Figures|4]and|5]show 
that H I and H2 emission originate in warm and hot regions respec- 
tively. The observed [N II] or [O II] emission cannot originate in 
the H2 region. Gas with a range of temperature and ionization must 
contribute to the observed emission. 

The Galactic ISIVI has two or more thermally distinct phases 
that are in pressure equilibrium (Field et al. 1969). Figure[T2lshows 
that gas with two temperatures but the same pressure exist for some 
values of the non-radiative heating rates. This occurs for cosmic ray 
fluxes ~ 4.2 dex times the Galactic background and for heating rate 
pre-coefficients of log ~ — 22.5. The range in gas properties for 
these two points is not broad enough to produce the observed range 
of emission however. For instance, in the cosmic-ray case the two 
hydrogen densities are ~ 10"' cm^-' and ~ 10^'^ cm^^, correspond- 

10^^ K and r ^ 10^^ K. Gas 



with this pair of densities, when co-added in the appropriate mix of 
filling factors, can account for some of the H I and H2 emission but 
none of the optical forbidden lines, as shown in Figure [TT] 

The full range of thermally stable clouds shown in Figure [T3] 
must exist and contribute to the spectrum. It is then necessary to 
adopt a weighting function, as described next, to decide how to 
co-add clouds with different conditions. This weighting function is 
not a physical model for why this mix of clouds exist. Rather it 
represents a parametric way of solving the spectroscopic problem 
by specifying how the volume of clouds depends on density. We 
will end up with an empirical fit to the spectrum and finally discuss 
various scenarios that might establish this range of conditions. 



5.2 Formalism 

We assume that the filaments are composed of cloudlets with a 
range of density n. The cumulative filling factor f{ii) of material 
with density below n (i.e. fractional volume filled by material with 
this density or less) is 
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where n is the total hydrogen density, a is a parameter to be de- 
termined by fitting the observed spectrum, and ^ is a normalization 
constant. The constant value 1 for densities above "high corresponds 
to the observation that the whole volume is filled by material of 
lower densities. The more commonly-used filli ng factor, described 
in Section 5.9 of lOsterbrock & Ferlandl |200^, is the dimension- 
less fraction of space that is filled with a particular phase of mate- 
rial. For a continuous density distribution, which we consider here, 
it is more useful to define a differential filling factor which will be 
a function of density. This is most easily defined in this cumulative 
form. 

The physical picture is that the aperture of the spectrometer 
takes in emitting material with a range of densities and correspond- 
ing temperatures but a single pressure. This might occur, for in- 
stance, if the non-radiative heating rate is not uniform in space or 
time. Gas would heat up or cool down, on the cooling timescale 
given in Figure[T4] in response to changes in the heating. The den- 
sity and temperature would change to maintain constant pressure 
but clouds with a variety of densities and temperature would con- 
tribute to the emission from a region on the sky. 

With these assumptions the line emissivity 'Xtzj (erg 
cm^-' s^') integrated over a distribution of clouds is given by 

\Ttj{n)df (15) 



Atzj 







dn 



(16) 



where the integration is over the physical limits to the filling factor, 
^ f{n) ^ 1, and the integration over filling factor is rewritten in 
terms of an integral over density. Here n\g^ and nhigh ars the lower 
and upper density bounds to the integration. In practice we simulate 
physical conditions and emission from clouds for a wide range of 
specific densities n and replace the integral with the sum 



/•=J 

;=0 / (=0 

where the density bins are separated by Sn/n -^OA dex. 



ing to kinetic temperatures of T : 



(17) 
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The hydrogen density averaged over the cloud distribution, n 
(cm^^), is similarly defined as 



n - 







ndf. 



This becomes 
n 



'^high "low 



(18) 



(19) 



With these definitions the luminosity of a line L/,„g ( erg s ' ) 
is related to the emissivity Atz j and the total volume containing gas 
V (cm^)by 



The mass contained in this volume is given by 



(20) 



(21) 



where fl is the mean mass per proton which we take to be lAu 
The line surface brightness is then 

12 



1.8704 X 10"'^ xATtjdl 



(22) 



The physical meaning of these quantities is as follows. The 
volume V is the total volume containing gas. The gas within this 
volume has a range of density and temperature and may be molec- 
ular, atomic, or ionized. The line emissivity An j takes into account 
both the atomic physics of the emission process and the fact that a 
line is only emitted for certain combinations of density and temper- 
ature. The weighting parameter a determines the amounts of gas at 
various densities. The density h depends on a since this sets how 
much weight is given to gas at various densities. It further depends 
on thermal stability criterion which sets nig„. 

There are three free parameters, the lower and upper density 
bounds and /thigh and the power-law index a. Physically mo- 
tivated limits are used to set the density bounds in Equation[T7] As 
described in Section |4!4l thermally unstable densities are excluded 
from the integrals. This is equivalent to assuming that gas that lies 
on unstable portions of the cooling curve, shown in Figure [T3] will 
quickly move to stable regions so that the volume is fully filled by 
stable gas. The CMB temperature is assumed to set a lower limit to 
the kinetic temperature which then sets the upper limit to the den- 
sity since the product nT is held fixed. This leaves only o; as a free 
parameter. 



5.3 Dependence on a 

The full emission-line spectrum, with many hundreds of thousands 
of lines, is computed for each point along the isobaric curve il- 
lustrated by Figure [T5] The emissivity of each line is stored and 
integrated over density using Equation [17] Integrated predictions 
are presented relative to Ha since it is easily observed and, of the 
bright optical lines, is the least affected by reddening. 

Figure[T5]shows that different lines have emissivities that peak 
in different regions along the isobaric line. Lines that are formed 
in warm or hot ionized gas, such as the [S II] or [O II] optical 
lines, will form in nearly the same region as Ha itself as shown 
in Figures [it] and [TS] We expect that the intensities of such lines 
relative to Ha will have little dependence on a but will depend on 
the chemical composition and the energy source. This is in contrast 
with molecular lines, such as the CO or H2 lines shown in Figure 
1151 which form in gas that is much colder and denser that the gas 
producing the H I emission. Because of the form of the filling factor 
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Figure 19. The integrated inten.sities of some strong lines are shown as 
a function of the power-law index a. The upper panel shows the H I Ha 
emissivity 47t j. The middle and lower panels show the predicted intensities 
of other lines relative to Ha for the cosmic ray (middle panel) and extra 
heating (lower panel) cases. 



function we expect that molecular lines will be stronger relative to 
Ha for a close to and that the intensity ratio decreases as a 
becomes more negative. 

Figure [19] shows the predicted line intensities as a function of 
a. Lines that form along with Ha, such as [S II] AA6716, are not 
strongly sensitive to the value of a. The line ratio is constant in 
the extra-heating case because Ha forms only in ionized gas along 
with the [S II] lines. The ratio decreases for large values of a in the 
cosmic-ray case because this gives strongest weight to the molecu- 
lar regions where the H I lines form by suprathermal excitation of 
H'' . For the regions where the line ratio is a constant the intensities 
reflect the abundances of the elements and the underlying heating 
process rather than the integration process. 

The intensities of molecular lines relative to Ha have a strong 
dependence on a. The power-law index a affects how the dense 
and cold molecular regions are added to the hotter ionized gas. Fig- 
urell9lshows that such lines as [O I] A63^.m, H2 2.12^i.m, and CO 
A863vtm have powerful dependencies on a. We will use these lines 
to determine a in the following section. 

Figure |20] shows the Ha emissivity and mean density as a 
function of a. The upper panel shows the emissivity 471 j of H I 
Ha as defined by Equation [T6] The emissivities are quite similar 
for when a is very small since this give greatest weight to the warm 
ionized regions where the emissivities are similar. They differ for 
large a since this weights the cold dense regions where suprather- 
mal excitation of H I is important in the cosmic-ray case. The total 
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Integration index a 

Figure 20. The upper panel shows the H I Ha emissivity An j as a function 
of the power-law index a. The lower panel shows the integrated density 
defined by Equation llSI 



volume of gas in all forms, molecule, atomic, and ionized, will be 
the observed luminosity of the Ha line divided by this emissivity, 
as given by Equation |20] The lower panel of Figure [20l shows the 
mean density defined by Equation[T8] The two cases differ slightly 
at low densities because the low-density cutoff given in Equation 
|19| is set by the thermal stability requirement. The total mass of gas 
in the filaments can be derived from the Ha luminosity and the in- 
formation given in this Figure together with Equation|2T] Note that 
although we give these results for Ha there is enough information 
in the Tables and Figures to do this with any emission line. 



5.4 The observed spectrum 

The next step is to use the observed intensities of molecular lines 
relative to Ha to set the power-law index a. This requires an ob- 
served spectrum that extends across the optical into the far-infrared. 
We take the Horseshoe region of the Perseus cluster as representa- 
tive. It is relatively isolated, at a projected distance of 25 kpc from 
the centre of NGC 1275, so confusion with other sources is min- 
imized. Its spectrum does not have underlying stellar absorption 
lines so optical H I lines can be measured without stellar contami- 
nation. 

The upper section of Table |5] gives the line ratios observed in 
the Horseshoe region together with predictions from the two mod- 
els. The optical (Gemini GMOS) and near IR (UKIRT CGS4) data 
are from Hatch et al ( 2006) and Hatch et al ( 2005) respectively and 
the mid-IR data are from ljohnstone et all ^20071) . All of the observed 
quantites in Table |5] have bee n corrected for an a ssumed Galactic 
extinction of E(B-V) = 0.163 dSchlegel et ailll998h . 

The lower section of Table |5] gives the observed Ha surface 
brightness, 5(Ha), and luminosity L(Ha) corresponding to the 
same region that the optical line ratios were extracted from. The 
averaged line emissivity, 47t j(lla), and averaged density, n, from 
the wave heated and cosmic-ray heated models are also listed. 



Table 5. Upper: Obse rved and predic t ed lin e strengths near the Horseshoe 
Position 11 region of IConselice et all i200ll) . Column 1 lists the emitting 
species, column 2 is the wavelength in microns, column 3 is the observed, 
extinction corrected flux ratio of each line with respect to Ha, column 4 
is the line flux ratio with respect to Ha predicted by the extra-heating 
model, column 5 is the line flux ratio with respect to Ha predicted by 
the cosmic-ray heating model. Note that the optical, NIR and MIR lines 
are each observed though a different aperture. See Figure [21] for a visu- 
alization of each aperture. Lower: Observed Ha surface brightness, S(Ha) 
( erg s^' cm^^ arcsec"^) and luminosity, L(Ha) ( erg s^' ) corresponding to 
the same region that the optical line ratios were extracted from. The inte- 
grated line emissivity averaged over the distribution of clouds, 4;ry(Ha), 
(ergs^' cm^') defined by equation 1 161 and the density averaged over the 
distribution of clouds, n, ( cm^') defined by equation lTsl in the wave heated 
and cosmic-ray heated models are also tabulated. 



Species 


A(|am) 


F(Une)/F(Ha) 


Pred./Obs. 


Pred./Obs. 






Observed 


Extra heating 


Cosmic ray 


a 






-2.10 


-0.35 


Hi 


0.4861 


0.24 


0.87 


1.04 


[Olll] 


0.5007 


<0.03 






[Nl] 


0.5199 


0.06 


1.88 


1.66 


He: 


0.5876 


0.04 


64 X 10-5 


0.37 


[0 I] 


0.6300 


0.19 


1.00 


0.93 


Hi 


0.6563 


1.00 


1.00 


1.00 


[Nil] 


0.6584 


0.71 


0.29 


0.34 


[Sll] 


0.6716 


0.23 


1.05 


1.40 


[Sll] 


0.6731 


0.17 


1.15 


1.08 


Hj 


1.957 


0.09 


1.28 


0.81 


Ht 


2.033 


0.03 


1.44 


1.41 


H2 


2.121 


0.08 


1.54 


1.02 


H2 


2.223 


0.02 


1.38 


0.77 


Hi 


12.28 


0.03 


1.00 


1.00 


[Ne II] 


12.81 


0.06 


4.6 X 10-3 


0.46 


[Ne III] 


15.55 


<0.02 






H2 


17.03 


0.06 


0.29 


0.53 


S(Ha) 




1 X 10-"* 






L(Ha) 




7 X 1035 






4;r/(Ha) 






9.40 X 10-21 


6.75 X 10-21 


n 






300 


2.63 X 10* 



The observational data presented in Table [5] were measured 
from a variety of ground and space based telescopes. These facili- 
ties each have different entra nce apertures, which we show overlaid 
on the Ha emission map of Conselice et a 3 l l200lh in Fig |21l As 
stressed bv l Johnstone et all ( 12007 ). these aperture corrections make 
it difficult to combine observations taken with different instrumen- 
tation. The 2^im NIR and optical spectra d Hatch etalll2005l.l2006l) 
were combined by scaled them to the same observed Ha flux. The 
Ha flux in the NIR data was estimated from the observed Pa line 
flux assuming that the Pa / Ha intensity ratio has the Case B value 
(8.45) for T = K. The MIR Spitzer observations of the low-y 
H2 lines were placed on the same flux scale as the optical/NIR lines 
by another method that did not assume Case B. Ha images of the 
region covered by the Spitzer entrance aperture were used to find 
the area covered by emitting gas and maps of the optical H I and 
MIR H2 emission wer e compared to introdu ce aperture corrections 
and derive flux ratios ( I Johnstone et al|[2007h . 

There are substantial uncertainties in the derived spectrum. 
The calculations presented above show that the Pa/Ha ratio is not 
expected to be given by Case B although the differences are not 
great (Table|4j. The extinction within the filaments is unknown so 
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Figure 21. Entrance apertures of the various instruments used to determine 
tlie line ratios presented in Table [5] and in the text. Gemini GMOS: long 
thin rectangle, UKIRT CGS4: small rectangle centred on brightest knot. 
The two larger overlapping rectangles are the positions of the Spitzer IRS 
SH aperture. Line fluxes were averaged between these two nod positions. 
The large circle is the beam of the IRAM (CO 1-0) observations. 



no correction is made. Renormalizing the optical / NIR H I spec- 
trum to agree with Case B intensities largely removes the effects of 
reddening. Given the uncertainties in both the observed and theo- 
retical predictions our goal is to reproduce the spectrum to within a 
fact or of two. 

ISalome et ^ ( l2008bl) detect the CO 1-0 line in the Horseshoe 
and note that the line is optically thick. This is in contrast to the 
optically thin optical and IR lines that are the focus of this paper. 
The observed flux of an optically thick line depends on the geome- 
try while that of an optically thin line depends only on the volume 
of e mitting gas. 

I Salome et 3 bOOSbh derive a mass of 10^ Mq for the re- 
gion of the Horseshoe included in their 1 1 arcsec radius beam (see 
Fig|21b. We can derive the total mass from the luminosity of any 
emission line using Equation |2T] We me asured the Ha lumin os- 
ity in the IRAM beam from the data of jConselice et alll200ll) as 
L(Ha) = 5.5x lO'^^ergs"'. 



5.5 The predicted spectrum 

We vary a to match the observed H2 12.28p.m/Ha intensity ratio of 
0.03 (Table|5j. This line ratio has several advantages. It has one of 
the widest ranges in excitation and ionization in the spectrum and 
so is a strong function of a. The IR H2 lines are, like the optical 
H I lines, optically thin so geometric effects do not enter. Finally, by 
comparing emission from an ion and molecule of the same element 
there are no additional uncertainties introduced by that fact that the 
composition is unknown. The uncertainty in the ratio is likely to be 
dominated by systematic errors which we cannot quantify. 

This procedure gives a = —0.35 for the cosmic ray and 
a — —2.10 for the extra heating cases. Both cases have negative a 
which, as Equation[T4]shows, means that the lowest-density gas has 
the largest volume. We are led to a picture in which the molecular 
gas lies in cold and dense cores while the low-ionization emission 
originates in extended low-density regions. The ionized material 
may be located in extended envelopes surrounding the molecular 
cores. 

The last two columns of Table [5] give the predicted line in- 
tensities relative to the observed values. A more complete set of 
predicted intensities is given in Appendix |A] Our goal is to match 
the spectrum to within a factor of two. Both heating cases do this 
for the majority of lines. The discriminants between the two cases 
are described in the next section. 

The last two columns of Table|4]give the integrated H I spectra 
for the two cases. The predicted Pa/Ha ratio is 0.0839 and 0.0992 
for the extra heating and cosmic ray cases, not greatly different 
from the Case B ratio of 0. 104. The Pa/Hj3 is 0.37 for both heating 
cases and is 0.29 for Case B. These differences are significant but 
the wide wavelength separation of the line pairs make them useful 
reddening indicators. 

5.6 The need for ionizing particles 

Several emission lines indicate that the cosmic-ray case best de- 
scribes the heating and ionization of the gas. Most lines are matched 
by both cases to within our desired factor of two. It is not possi- 
ble to formally quantify a goodness of fit in Table |5] because the 
uncertainties are dominated by systematic errors which we cannot 
quantify. Rather we point to a few lines which clearly indicate that 
the cosmic-ray case applies. 

Figures [17] and [TS] show the critical distinction between the 
two cases. The extra-heating case behaves as expected for a 
coUisionally-excited thermal gas. As the temperature increases, go- 
ing from left to right in the figures, we see lines of [O I], [N I], 
[Nil], [On], [Olll], He I, and eventually [Nelll], increasing in 
intensity. This corresponds to increasing ionization and excitation 
potentials as expected for a collisionally-ionized gas. As stressed in 
the discussion of Figure[8]the distribution of ions and molecules in 
the extra-heating case is sharply defined by temperature. At a given 
temperature hydrogen will be nearly all ionized, atomic, or molec- 
ular with the change from one species to the next occurring over a 
relatively small range in temperature. 

This behavior is to be contrasted with the cosmic-ray case. 
Here a certain level of ionization and excitation is present at all gas 
kinetic temperatures because the cosmic-ray ionization rates are 
only indirectly related to temperature. The result is that some H" 
is present in regions that are predominantly H2 or H+. We find the 
somewhat paradoxical result that the higher-ionization (3.00 Ryd 
to create) [Ne III] A3869A line becomes strong at lower tempera- 
tures than the lower-ionization (2.58 Ryd) [O III] A15007,4959A 
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lines. This is because of the importance of charge-transfer recom- 
bination in setting the ionization distribution of O. The process 
0+°-hH° ^ 0+""' +H+ is very fast for both 0+ and 0+^ but 
is slow for He+, Ne+, and Ne+^. Larger amounts of are present 
in the cosmic-ray case because of the mixing of different ioniza- 
tion species in the same region. The result is that the cosmic-ray 
case produces significant [Ne II] Al2.81ii-m and [Ne III] A3869A 
emission without exceeding the upper limit to [O III] A5007A. 

He I /15876A provides another distinction between the cos- 
mic ray and extra-heat cases. He" is more efficiently ionized by 
suprathermal than by lower-energy thermal particles due to its large 
ionization potential. The A5876A line is in fair agreement with the 
observed value in the cosmic ray case while it is underpredicted by 
more than 4 dex in the extra heating case. 

Figure [TS] shows the fundamental tension experienced by the 
two cases. We do not include the thermally unstable temperatures 
indicated by the grey bar in the Figure. The extra-heating case 
could only produce significant He I and [Ne III] emission by in- 
cluding large amounts of unstable gas. This would overpredict the 
unobserved [O III] A5007A line by large factors. The cosmic-ray 
case produces strong He and Ne emission while not overproducing 
[OiiiL 

Ijohnstone et all ( |2007|) discuss the puzzle posed by the ob- 
served [Ne III] emission from filaments that do not show signifi- 
cant [O III] X5007A. Neon is special because of the high ionization 
potentials compared with other second row elements like N or O. 
[Ne III] was not detected in the Horseshoe but was seen in another 
region where [O III] was weak. Ijohnstone et all ( l2007h suggested 
that a low kinetic temperature in the ionized gas, perhaps the result 
of enhanced abundances and cooling efficiency, might be the cause. 
Cosmic-ray ionization is another possibility. 

These calculations demonstrate that cosmic rays are needed 
to account for the spectrum. Suprathermal particles such as those 
produced by cosmic rays create a gas with a broad range of ion- 
ization at any one point. Charge transfer recombination then pro- 
duces the observed unusual ionization distribution of the heavy ele- 
ments. Although we have explicitly considered cosmic rays, which 
are known to exist in the cluster environment, any source of ion- 
izing particles would work as well. Although gentle heating by 
MHD waves is ruled out, MHD waves could still be the funda- 
mental source of the ionization since high-energy particles may be 
prod uced in situ by MHD phenomena such as magnetic reconnec- 
tion /Lazarian''2005'). 

Crawford & Fa bianI j 19921) showed that photoionization by a 
very hard continuum source can produce the optical spectrum. 
High-energy photons produce primary and Auger electrons that 
have very high energies and so behave as ionizing particles. Their 
model worked for the same reasons that the cosmic-ray model dis- 
cussed above works. The ionizing particles generate a gas with a 
mix of molecules, atoms and ions, charge transfer effects are impor- 
tant, and a peculiar low-ionization spectrum is produced. Photoion- 
ization by very hot objects is another viable source of the deduced 
population of suprathermal particles. 

A form of mixing layers might also produce the ionizing par- 
ticles. If hot ionized gas were to mix with cold molecular gas the 
resulting physical state would be determined by the microphysics 
of the coUisional interactions rather than by the ratio of kinetic tem- 
peratures. Any energy source that deposits its energy as ionizing 
particles rather than as heat can produce the observed spectrum. 
Our predictions do not depend on the fundamental energy injec- 
tion, only on the eventual production of a population of suprather- 
mal particles. We shall refer to this as ionizing particle heating in 



the remainder of this paper to reflect the fact that there are many 
possible sources of these particles. 



6 DISCUSSION 
6.1 The geometry 

Our deduced value of a = —0.35 corresponds to the majority of the 
volume being filled by low-density ionized gas. The density of the 
molecular regions is typically n„oig/nig„ ~ lO"* times higher than 
the density of the ionized gas (Figure [Tst. The ratio of volumes is 
then Vmole/^ion 10"'" ~ 0.09. By contrast, most of the mass is in 
dense molecular cores. The mass depends on w"+^ so the ratio of 
masses is M„joig / Mion ~ 10"''"+'^ ~ 90. The small dense molecular 
cores may provide the reservoir of material for an extended halo of 
low-density ionized gas. If the gas is constrained to follow magnetic 
field lines then the individual clouds may resemble comets with a 
cold dense core and extended warm low-density tail. 

The observed surface brightness in Ha and the predicted 
4;rj(Ho!), both listed in TablejS] can be combined with Equation 
l22lto find the line-of-sight thickness of emitting material. The de- 
duced thickness is rfZ = 1.0 x lO"^ cm = 0.3 pc. Note that this is the 
thickness of the entire cloud, not just the H+ region, although most 
of this is filled by ionized gas. 

This depth is similar to that found bv lFabian et al( l2003h from 
fundamentally different assumptions. They assumed that the opti- 
cal H I lines are formed by recombination with Case B intensities, 
that the pressure of the H I-emitting gas was equal to that of the sur- 
rounding hot gas, and that the gas was fully ionized with a temper- 
ature of ~ lO'^K. Ifthe gas is excited by ionizing particles then H I 
is far more emissive than predicted by Case B due to the enhance- 
ment by coUisional excitation (Table |4j. This large enhancement 
makes up for the fact that hydrogen is not fully ionized across the 
volume so the H I emission measure, the product iierip dl, is much 
smaller. These two factors nearly cancel leading to a similar depth. 

For this surface brightness argument to work the entire en- 
trance aperture of the spectrometer must be covered by emitting 
gas. A slit width of 1 " corresponds to a projected width of ~ 350 pc 
at the distance of Perseus. Each cloud is --^ 0.3 pc thick along the 
line of sight but must cover a region of the sky --^ 350 pc across. The 
aspect ratio of the width to depth of the emitting gas is > lO' : 1 
since the region of the Horseshoe we observed is well resolved on 
the sky. This aspect ratio, although surprising, is similar to other 
geometric quantities in the system. HST observations show that the 
projected width of the narrowest parts of the filaments is ~ 70 pc 
jFabian et alll2008h while they have lengths of tens of kpc. The fil- 
aments are composed of remarkably thin and long clouds. 

The filaments are most likely composed of a network of much 
smaller threads. In the magnetic-field dominated constant-pressure 
case assumed here gas is free to move along but not across field 
lines. The line surface brightness suggests that individual threads 
have a width of a fraction of a pc. We view this geometry from a 
direction that is nearly orthogonal to the field lines. This is an ob- 
servational selection effect introduced by the way that isolated fila- 
ments were chosen for study. The narrowest filaments, with a width 
of 70 pc, must be composed of --^ 10^ — lO-' of these thin threads. 
If the threads are uniformly distributed across the filaments, and 
the filaments have a depth equal to their observed width, then the 
mean separation between threads is ~ 10 pc, about thirty times their 
thickness. 

There are two possible physical interpretations for our a 
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power-law index. It could represent a radial gradient in density 
across a thread normal to the field lines. In this case the dense 
molecular regions would form a very small linear core with the 
lower-density ionized gas surrounding it in a much larger cylin- 
der. We consider this unlikely since gas could not move across field 
lines to change its density in response to changes in the heating. The 
more likely possibility is that the a index represents the variation of 
the density along magnetic field lines. In this picture the ionizing- 
particle density would vary with position or time, the temperature 
of the gas would adjust itself in response, and then expand or con- 
tract to maintain constant pressure. A magnetic field line would 
then have most of its volume occupied by lower-density ionized gas 
intermixed with dense molecular regions. The cold dense molecu- 
lar regions could become warm and ionized, and return to being 
cold and molecular, as the heating rate changes. 

6.2 Internal extinction across a thread 

The hydrogen column density through a thread normal to its surface 
is A'(H) iidl ~ 10^^'^ cm^^ . If the gas has a Galactic dust to gas 
ratio, as we have assumed, then the total visual extinction across a 
thread will be A,, --^ 17 mag. The effects of this extinction depends 
on the geometry. The dense molecular regions, having most of the 
column density, would also be most extinguished. This region de- 
tected emits MIR and sub-mm emission where the extinction will 
have little effect. The column density through the partially ionized 
gas which produces most of the optical emission is small so this 
emission is not directly affected by internal extinction. We could 
not detect optical emission from ionized regions which happened 
to lie on the far side of the molecular cores of course. 

6.3 The total mass at position 11 

The total mass can be derived from the luminosity of any emission 
line using the formalism derived in Section [5^ Section [5!4l derived 
an L(Hot) of 5.5 x 10*^^ erg s^' for position 1 1 and the beam used 
bv lSalomeet"3 ( l2008bl) . Equation 12 II gives a total mass of 2.1 x 
10*^ Mq. This is in reasonable a greement with the mass of 10*^ Ma , 
derived bv lSalome et ail ( l2008bl) from observations of the CO 1-0 
transition. 

The mass of 2. 1 x 10^ Mq derived from emission lines counts 
only the spectroscopically active gas. As Figure [TS] shows, little 
emission is produced by cold gas at the high-density end of the dis- 
tribution. Gas denser than our high-density limit cannot exist and 
be in pressure equilibrium because the CMB sets the lowest possi- 
ble temperature. However there is no observational constraint that 
actually limits the amount of material that could be at this highest 
density and the CMB temperature. A large reservoir of very cold 
high density gas could exist and remain undetected. Our derived 
mass is a lower limit as a result. Gas at the very coldest limit of 
the distribution could serve as the source of fresh material for the 
spectroscopically active regions. 

6.4 The pressure of the molecular cores 

A number of studies have estimated the pressure of the molecular 
gas. Next we compare these with our derived temperature and den- 
sity which assume a pressure equal to that of the surrounding hot 
gas. 

Ijaffe et all ( 1200 l l) found a density )Ih > 10^ cm"^ and T ~ 
10^'^ K from the distribution of H2 level populations. The gas pres- 



sure, nT > W^-^ K, was substantially higher than the pressure of 
the surrounding hot gas. For our best fitting ionizing particle model 
most of the v = 0, low-7, H2 emission originates in gas with a tem- 
perature of about T ~ 10' "' K and a density ;jh ~ 10^'^ cm^^, as 
shown in Figure[T5] Excited-v emission occurs at a higher temper- 
ature, T --^ 10^'^ K, and lower density, njj ~ lO-''^ cm^-'. We do not 
present level-excitation diagrams in this paper but note that we have 
fitted the observed H2 spectrum, from which the H2 level-excitation 
populations were derived, to within a factor of two. Our model pro- 
duces a population distribution that is consistent with these obser- 
vations. 

ISalomeetall l l2008iJ) argue that the filaments are composed of 
an ensemble of Galactic Giant Molecular Cloud-like components. 
By analogy with GMCs they argue that the density in the CO re- 
gion is several hundred particles per cubic centimeter and that the 
gas has a kinetic temperature of T ~ 20 — 40 K. The resulting pres- 
sure, riT ~ lO"' "^ cm^-' K, is 3.2 dex smaller than the pressure of the 
surrounding hot gas. 

We avoid geometrical complexities by focusing on optically 
thin emission lines. Figure [Tsl shows that the optically thin emis- 
sivity of the CO A 863 [im line peaks at Tj^qq ~ 20K and n ~ 
10^^ cm~^ . This tempe r ature i s similar to the CO temperature 
quoted by ISalome et 

M ( l2008al) but the density is considerably 
higher in keeping with our assumed high gas pressure. The dif- 
ference is due to the higher cosmic ray density in our simulations. 
Well-shielded regions of Galactic GMCs are similarly heated by 
cosmic rays but with a considerably lower intensity equal to the 
Galactic background. The cosmic ray density that produced the 
peak in our simulations is ~ 10* times the Galactic background. 
This accounts for the GMC-like temperatures at such a high pres- 
sure. It is not possible to make a specific prediction of the CO 1-0 
brightness temperature without first deriving an explicit geometric 
model of the cloud, a task beyond the scope of this paper. 

We conclude that both the H2 spectrum and CO bright- 
ness temperature are consistent with a gas pressure of nT = 
cm ^ and excitation by ionizing particles. 

6.5 Gravitational stability 

The stability of the material in the filaments to gravitational col- 
lapse may be influenced by a number of factors: basic geometry, 
turbulence and magnetic fields. 

If the support of the filaments is purely by thermal pressure, 
then they will be liable to collapse when their size becomes greater 
than the Jeans length, 

~ lO^^cgnj '/^cm, (23) 

where C5 is the sound speed (or effective virial velocity) in kms^' 
and the density is lO^nscm^-'. For cg = 1, 05 = 1 the thickness of 
the molecular material will be 0.3 pc, a similar scale to that inferred 
for the Ha emitting region. We infer above that the molecular mate- 
rial has a volume filling fraction about one tenth that of the ionized 
gas, so the size scale of the molecular material will be comparable 
to this. This suggests that the gravitational stability of the molecular 
component may have a significant role in determining the structure 
of the filaments. 

Magnetic fields and turbulence may play a significant role in 
supporting the molecular material against gravitational collapse. 
The standard criterion for support against collapse by uniform mag- 
netic fields is to compare the column density along a flux line to 
the critical value B/2kG^^^. If the dense material is threaded along 
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flux tubes, this is not a significant constraint as is stands. Likely 
more important is tiie coupled effects of turbulence and randomly 
directed internal magnetic fields. These will act to increase the ef- 
fective virial velocity of the material, and so support it against col- 
lapse. 

However, both turbulence and random magnetic fields are in- 
herently transient, and so must be driven in some fashion to main- 
tain their role in supporting the material. In the local interstellar 
medium, this driving is achieved by a number of processes, in- 
cluding star formation. While there is some evidence of star for- 
mation in the filaments, this appears to be a relatively localised 
phenomenon. 

What will happen if part of the molecular material in a fila- 
ment starts to undergo a gravitational collapse? We find that the 
cooling timescale (~ 10^ yr) is significantly smaller than the dy- 
namical timescale (~ 3 x 10^ yr), so the material may initially frag- 
ment. However, the enhanced density will mean that the net tidal 
force between the material and the surrounding diffuse envelope 
will increase, both allowing the filament to break up, and enhanc- 
ing the net heating. This process may thus naturally maintain the 
filaments with marginal gravitational stability, in a similar fashion 
to that in which galactic disks are maintained with a Toomre Q pa- 
rameter around unity. Note that if this is the case, it requires that the 
heating rate is determined by local processes, such as drag-driven 
reconnection. 



6.6 The energy budget 

The energy dissipation timescale for the line-emitting material is 
short, requiring that the energy be replenished in situ. While it is 
clear from Figure [TSl that below lO'^K the detailed cooling mech- 
anisms depend on the temperature, from Figure [T3] the overall 
cooling rates do not depend strongly on these details. The heat- 
ing need to power the emission regions in steady state ranges from 
A ~ 10^^^*' ergcm^-' s^' while the thermal timescale of the gas 
is ~ 10^ yr (Figure[T4]l. 

A lower limit to the energy density of MHD waves in equipar- 
tition with a magnetic field of 100 jlG, the weakest B capable of 
confining the clouds, is 



= 4.4x 10"'"(S/100mG)^ ergcm" 



For the quoted heating rate this energy will be dissipated in a 
timescale 

We = t/wave/A= 10^*' (S/lOO/iG)^ yr 

over which the wave energy has to be replenished. 

The timescales for the cosmic-ray heated case are somewhat 
longer. The Galactic background cosmic -ray density corresponds 
to an energy density of ~ l.SeVcm^-' (Webber 1998), mostly in 
high-energy nuclei. The cosmic-ray intensity required to maintain 
the material at nT --^ lO*'^ cm^-'K ranges between 10^-10* times 
the Galactic background value, depending on the filament density 
(Figure [T2}. Molecular regions are predicted to have a density of 
~ 10^ cm^-' and require a cosmic-ray --^ lO-* times the back- 
ground. The low-ionization gas requires a rate ~ 10* times the 
background with a density of «h ~ 10^'^ cm^-*. For these values 
the energy dissipation time for the molecular and low-ionization 
regions are ^ 10"' yr and ~ 5 x 10^ yr respectively. 

These replenishment timescales are generally of order or 
longer than the cooling time (Figure [74l (, the time required for the 
gas to adjust itself to a changing heating rate. This means that non- 
equilibrium effects, in which the ionization, heating, and cooling 



decouple from one another, may be important. The short heating 
/ cooling times suggests that the gas pressure can fluctuate about 
its equilibrium value in response to variations in the environmental 
heating, and may drive sound waves. The resulting structures may 
be intrinsically quite dynamic and disequilibrium. 

The size scale of a thread appears to be ~ 0.3 pc. Assuming 
that they are cylindrical, then if they are heated through their sur- 
faces the energy flux required to keep them in equilibrium is 

'■fil 



-A~0.05ergcm"2s"'. 



(24) 



If the field which provides this energy has an energy density X 
times the thermal energy density inferred for the surrounding X-ray 
bright gas, the net speed of energy transport through the surface of 
the filaments must be 
'■fil -1 



10° 



(25) 



A'0.3pc 

While this suggests that it is easier to heat more narrow filaments, 
the same energetic constraints will apply to an assembly of fila- 
ments if the observed surface brightness is to be produced. For 
comparison, the Alfven velocity is 



B 



7x \(fn^^/^ (B/lOO/xG) cms"', 



(26) 
(27) 



and the sound speed will also be of order Ikms"'. This suggests 
that the energy required to heat the filaments can only be trans- 
ported through the molecular material by non-thermal particles 
or radiation. However, this does not rule out it being converted 
into this form locally to the filaments, as energy transport is eas- 
ier through the more diffuse components with higher characteristic 
transport speeds. 

The filament material does not appear to be supported against 
gravity by lateral kinetic energy, as the observed speeds are below 
the local Keplerian velocity, and the filaments are being dragged 
away from the stellar populations which are believed to have 
formed within them. If this is the case, the net rate of drag heat- 
ing of the material per unit mass will be approximately 



- tidal 



I'Kep^'drag 



0.3ergg 's ' 



:5x 10"^"n5ergcm"^s"' 



(28) 
(29) 



assuming VKep — 10 kms 



' 5 ^'drag — 1 0^ km s ' and that the mate- 
rial has a density of lO^nscm"'. This energy input rate is in the 
range of interest for heating the filaments, as alre ady suggested 
from the overall energy budget bv lPope et all ( l2008h . Note that the 
spatial distribution of this energy input will be determined by the 
details of the drag processes between the phases, so that the result- 
ing energy input may appear either in the dense material or within 
the halo of the material. 

Note that these scaling values can be formed into a simple 
estimate of the rate at which mass passes through the system of 
filaments, by dividing the total mass inferred for the molecular fil- 
aments by the dynamical timescale of the material, i.e. 

=^McoVd,ag/'-^50MQyr-l, (30) 

which is similar to the mass processing rate inferred for the radia- 
tive cooling of the hot material in these regions. 

At high Reynolds number, the drag force is roughly equal to 
the ram pressure between the phases. At the relative density of the 
X-ray hot medium and the molecular gas, this would suggest that to 
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prevent the dense material moving at free fall velocities, the relative 
motion of the hot material must be at (imol/"hot)'^^i'ff> which are 
not likely to be the case. The filament material may be transitory 
on a timescale ^ vgi/g ~ 3 x 10^ yr, in which time the filament 
material will travel --^ 1 /2(v/vK:ep)^ of the distance to the core of the 
cluster. Alternatively, their effective cross section may be increased 
either by an effective viscosity or by magnetic fields. 

The interphase drag will lead to magnetic field lines in the 
diffuse phases becoming wrapped around the dense material. This 
suggests that the magnetic structure may be somewhat like the bub- 
ble in the solar wind maintained by the earth's magnetosphere. 
The resulting magnetic segregation would limit high energy par- 
ticle transport into the molecular material. It is likely that some 
high energy particles will be able to leak in, either through mag- 
netic neutral points, as a result of turbulent motions, or by diffusion 
resulting from field inhomogeneities. Alternatively, the drag-driven 
turbulence may itself be the source of the high-energy particle flux, 
as a result of reconnection events. 

Overall, it is clear that the dissipation of the energy resulting 
from buoyancy has potential as a source of heating for the fila- 
ments. However, if this is the source of the heating, the means by 
which the energy is transformed into the non-thermal ionizations 
which we infer are the immediate means by which the molecular 
material is heated must be complex. 

6.7 What is the source of the ionizing particles? 

The chemical-ionization state of the gas strongly suggests that the 
the heating and ionization is produced by ionizing particles. Several 
considerations suggest that the population of ionizing particles is 
not produced by a strongly enhanced version of Galactic cosmic 
rays. 

The energy density Galactic cosmic rays is dominated by 
particles with energies around IGeV in the solar neighbor- 
hood (where at higher energies the energy spectrum is oc E^^'^, 
[Spitzer & Tomasko ( 1968)). Most of the energy is in high-energy 
protons which have a long range and do not interact strongly with 
matter. The CR particles lose energy by ionizing the material they 
pass through, primarily by producing secondary electrons with en- 
ergies around 50eV. These secondary electrons further ionize and 
heat the gas, but they will be transported by relatively small dis- 
tances. 

As a rough estimate of their range, the vel ocity of CR protons 
at 1 Ge V is V = a/3/2c, so from equation (3) of lSpitzer & Tomaskd 
( Il968 l) the interaction cross section with atomic hydrogen is ~ 
10^'^ cm^. Each interaction causes the particle to lose ~ 50eV of 
its energy. The range of the CR heating particles is therefore 

LcR=^— ^^700«5-'pc. (31) 
noAE ^ 

However, this range may be significantly reduced by the magnetic 
field. The gyroradius of the protons in a 100 /iG = lOnT field is 

rg = ^ = 3x10"' 7)3 cm (32) 

so the protons will be tightly coupled to the magnetic field lines, but 
(at least at first order) free to travel along them. Non-uniformities 
in the field geometry will allow particles to diffuse across the field 
lines. 

The distance which the particles can travel in a dissipation 
time is ~ 30pc in the molecular gas and ~ 150kpc in the atomic 
gas. This suggests that the cosmic rays can travel with little reduc- 
tion in intensity across the size scale of an individual thread. It is 



difficult to see how the flux of cosmic rays could be modulated to 
the extent that we need to produce the range in density shown in 
Figure[T2l 

Even if we were to assume that the large flux of CR required 
to heat the ionized material could be excluded from the molecu- 
lar core, this would result on a pressure acting on the molecular 
material which would be inconsistent with our models. The Galac- 
tic background cosmic ray energy density of --^ l.SeVcm^^ cor- 
responds to a pressure of ~ 2 x 10"* cm^^ K. Our spectral model 
requires a range of ionization rates between 10^ — 10^ times 
the background. The corresponding pressure is between lO'' -' — 
lO'*' "' cm^-' K, with the top of this range being far larger than that 
inferred for the molecular material. The logical conclusion is that 
the inferred suprathermal ionization is not associated with a popu- 
lation of relativistic nuclei that is similar to Galactic cosmic rays. 

As described above, other possible sources of suprathermal 
electrons include magnetic reconnection, mixing layers, and pho- 
toionization by very hard radiation fields. 



6.8 The history of the gas 

The origin of the filaments is unknown. Two possibilities are that 
molecular gas was ejected from the ISM of the central galaxy or 
that they formed f rom the surrounding hot gas which fills the cluster 
( iRevaz et alll2008l) . 

We have shown that the filaments are energized by a popula- 
tion of ionizing particles. High-energy particles destroy molecules 
in addition to heating the gas. Figure [15] shows the the peak CO 
and H2 line emission occurs at « ~ lO^ -' cm^-' while Figure [T2I 
shows that this occurs at a cosmic r ay ionization rate of ~ 10^'^ 
times the Galactic background, which lwilliams et ail( ll998h give as 
5 X 10^'^ s^'. The corresponding H2 dissociation rate of ~ 1.6 x 
lO'"* s^' corresponds to an H2 survival time of only --^ 2 x 10^ yr. 
This is very close to the H2 formation timescale for the dusty case, 
as give in Section 13.51 above, but far shorter than the formation 
timescale in the dust-free case. We conclude that that gas must be 
dusty to support the observed molecular inventory against rapid de- 
struction by ionizing particles. 

Perhaps surprisingly, the abundances are fairly close to that 
of the ISM of our galaxy. In photoionization equilibrium the 'ther- 
mostat effect' governs the strength of forbidden lines relative to 
hydrogen lines. Their intensit y ratio has no direct dependence on 
the abundances ( lFerlanj|2003h . In the non-radiatively heated cases 
considered here the strength of the forbidden lines is a product of 
both the temperature of the gas and the abundance of the element. 
In this case the intensities of forbidden lines relative to hydrogen 
lines depends linearly on the abundance of the heavy element rela- 
tive to hydrogen. This suggests that the metallicity of the filaments 
is within a factor of two of solar. 

Th e [Ca II] doublet n ear A7320A is strong if grain are not 
present (Kingdon et al 1995). The absence of this line shows that 
Ca is depleted by at least an order of magnitude, also suggesting 
that dust is present. 



7 CONCLUSIONS 

(1) We propose a model in which the filaments have a constant 
gas pressure that is set by the pressure of the surrounding hot gas. 
The overall geometry is set by magnetic field lines. Gas is freed to 
move along field lines but not across them. 
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(2) We consider gas heating by dissipation of MHD wave en- 
ergy and by cosmic rays. If a range of fieating rates occur then the 
gas will expand or contract along field lines to maintain constant 
pressure. Low heating rates result in dense molecular cores while 
high heating rates produce an extended warm ionized gas. 

(3) The relative proportions of molecular and low-ionization 
emission is set by the proportion of high and low density gas that is 
present. In our final best-fitting model most of the volume is filled 
by low-density ionized gas while most of the mass is in dense cold 
molecular cores. 

(4) Both wave and cosmic -ray heating can match most emis- 
sion line intensities within a factor of two. The distribution of low- 
ionization emission lines shows that charge transfer strongly affects 
the ionization. This is only produced in the cosmic-ray case, where 
a population of suprathermal electrons produces a gas with a wide 
mix of ionization. Although we explicitly consider cosmic rays any 
source of ionizing particles would work as well. 

(5) The emission measure shows that the filaments are com- 
posed of a network of much smaller threads, each 0.3 pc in ra- 
dius and separated by --^ 10 pc. A single filament is composed of 
hundreds to thousands of these threads. 

(6) The mass we derive from fitting the entire spectrum is in 
reasonable agreement with mass estimates from optically thick CO 
lines. There may be an even larger reservoir of dense molecular 
gas at the CMB temperature which would not be spectroscopically 
active. 

(7) The observed molecular inventory requires that dust be 
present to sustain the needed H2 formation rate. The presence of 
dust suggests that the filaments have not been strongly shocked and 
did not form out of the surrounding hot and presumably dust-free 
gas. 

(8) The metallicity is within a factor of two of solar. Both this 
and the presence of dust suggests that the gas originated in the ISM 
of the central galaxy. 

(9) The calculations presented here use the development version 
of Cloudy, which we have expanded to include new H2 collision 
rates and resolve H" configurations into nl terms. The critical densi- 
ties of the H2 levels which produce the observed emission are con- 
siderably lower than previously estimated. The ionizing particle- 
excited H I spectrum has relative intensities that are not too differ- 
ent from Case B but has emissivities which are far higher. 
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Table Al. Predicted UV and optical line ratio.s. 



Line A { yun) Extra Heat Cosmic Ray 



Mg II 


0.2798 


1.64E-01 


1.04E-01 


[On] 


0.3727 


1.18E+00 


1.27E+00 


Hi 


0.3798 


6.41E-03 


1.09E-02 


Hi 


0.3835 


9.38E-03 


1.52E-02 


[Ne III] 


0.3869 


1.28E-10 


1.32E-02 


He I 


0.3889 


7.45E-06 


9.01E-03 


HI 


0.3889 


1.45E-02 


2.24E-02 


Hi 


0.3970 


2.44E-02 


3.52E-02 


[S II] 


0.4074 


4.78E-02 


4.16E-02 


Hi 


0.4102 


4.59E-02 


6.08E-02 


Hi 


0.4340 


1.02E-01 


1.20E-01 


He I 


0.4471 


8.51E-07 


5.39E-03 


Hell 


0.4686 


1.40E-13 


4.36E-03 


Hi 


0.4861 


2.27E-01 


2.71E-01 


vFv 


0.4885 


2.39E+00 


1.43E+00 


[0 III] 


0.5007 


9.73E-10 


2.01E-03 


[Nl] 


0.5199 


1.13E-01 


1.02E-01 


[Nil] 


0.5755 


7.15E-03 


7.24E-03 


He I 


0.5876 


3.19E-06 


1.83E-02 


[O I] 


0.6300 


1.91E-01 


1.76E-01 


Hi 


0.6563 


l.OOE+00 


l.OOE+00 


[Nil] 


0.6584 


2.15E-01 


2.49E-01 


He I 


0.6678 


1.84E-06 


5.41E-03 


[Sll] 


0.6716 


2.30E-01 


2.37E-01 


[Sii] 


0.6731 


1.72E-01 


1.78E-01 


He I 


0.7065 


8.36E-06 


1.75E-03 


[At III] 


0.7135 


1.77E-06 


8.57E-03 


[Ca II] 


0.7306 


9.39E-04 


2.17E-03 


[Oil] 


0.7325 


4.48E-02 


4.95E-02 


[S III] 


0.9532 


4.90E-05 


8.43E-03 
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APPENDIX A: THE EMISSION-LINE SPECTRUM 

Table lAll lists the intensities of some of the brighter optical emis- 
sion lines while Tables IA2I and IA3I give the predicted lines in the 
near IR through sub-mm regions. 
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Table A2. Predicted near IR line ratios. 



Table A3. Predicted mid IR - sub-mm line ratios. 



Line 


A(|am) 


Extra Heat 


Cosmic Ray 


He I 


1.083 


7.72E-05 


2.37E-02 




1.631 


4.45E-06 


2.64E-04 


H2 


1.636 


5.55E-07 


1.63E-04 




1.639 


9.00E-06 


2.00E-04 


H-> 


1.645 


1.62E-06 


2.58E-04 


Ht 
^^z 


1.650 


2.05E-04 


4.82E-03 


^*Z 


1.658 


1.83E-07 


2.29E-04 


H-> 
^*z 


1.666 


2.80E-04 


4.77E-03 


^*z 


1.675 


6.16E-07 


1.61E-03 


^^z 


1.687 


2.58E-03 


3.13E-02 


^*z 


1.695 


2.28E-07 


5.78E-04 


^^z 


1.714 


2.32E-03 


1.68E-02 


^^z 


1.719 


6.73E-07 


6.18E-04 


^^z 


1.747 


2.43E-07 


3.69E-04 


^^z 


1.748 


1.58E-02 


6.32E-02 




1.748 


1 58F-02 


6 37F-02 


H-> 
^*z 


1.780 


5.95E-07 


3.60E-04 


^^z 


1.787 


1.05E-02 


4.35E-02 


^^z 


1.835 


5.57E-02 


6.37E-02 


H I 


1.875 


8.39E-02 


9.92E-02 


^*z 


1.891 


2.83E-02 


6.50E-02 


^*z 


1.957 


1.15E-01 


7.26E-02 


H-> 
^*z 


2.033 


4.33E-02 


4.24E-02 


^*z 


2.121 


1.23E-01 


8.17E-02 


H I 


2.166 


4.84E-03 


6.35E-03 


^*z 


2.201 


1.50E-03 


3.16E-03 


^^z 


2.223 


2.76E-02 


6.17E-02 


^^z 


2.247 


1.12E-02 


2.46E-02 


^*z 


2.286 


5.77E-04 


6.93E-03 


H-> 
^*z 


2.406 


9.16E-02 


7.88E-02 


^*z 


2.413 


3.05E-02 


6.81E-02 




2.423 




5.73E-02 


^*z 


2.437 


2.41E-02 


2.35E-02 


^*z 


2.454 


5.57E-02 


3.50E-02 


^^z 


2.475 


1.26E-02 


2.90E-02 


^^z 


2.499 


2.42E-02 


2.77E-02 


H2 


2.527 


4.69E-03 


1.94E-02 


H2 


2.550 


8.45E-03 


2.52E-02 


H2 


2.558 


3.09E-03 


9.46E-03 


H2 


2.595 


1.30E-03 


9.44E-03 


H2 


2.626 


2.18E-02 


4.93E-02 


H2 


2.710 


9.63E-04 


2.06E-02 


H2 


2.718 


3.59E-04 


5.05E-03 


H2 


2.719 


1.27E-03 


8.02E-03 


H2 


2.730 


1.03E-03 


7.14E-03 


H2 


2.802 


7.75E-02 


6.66E-02 



Line A ( |am) Extra Heat Cosmic Ray 



H2 


5.051 


1.66E-02 


3.70E-02 


H2 


5.510 


1.06E-01 


1.57E-01 


H2 


5.810 


2.89E-03 


1.16E-02 


H2 


6.107 


6.63E-02 


7.30E-02 


H2 


6.907 


2.97E-01 


2.47E-01 


H2 


8.024 


1.04E-01 


1.04E-01 


H2 


9.662 


2.16E-01 


1.83E-01 


H2 


12.28 


2.97E-02 


3.15E-02 


[Ne II] 


12.81 


2.74E-04 


2.78E-02 


[Ne III] 


15.55 


3.17E-10 


8.17E-03 


H2 


17.03 


1.79E-02 


1.94E-02 


[S III] 


18.67 


1.48E-05 


3.72E-03 


H2 


28.21 


3.12E-04 


2.07E-03 


[S III] 


33.47 


2.89E-05 


7.82E-03 


[Sim] 


34.81 


2.30E-02 


8.98E-01 


[Ne III] 


36.01 


2.76E-11 


7.31E-04 


[O III] 


51.80 


8.15E-11 


1.29E-04 


[01] 


63.17 


7.52E-02 


2.75E+00 


[Sil] 


68.40 


1.43E-04 


4.06E-03 


[Oiii] 


88.33 


1.16E-10 


1.81E-04 


[Nil] 


121.7 


4.79E-03 


7.83E-03 


[01] 


145.5 


7.48E-03 


1.43E-01 


[Cll] 


157.6 


2.57E-02 


1.31E-01 


CO 


323.6 


4.17E-04 


4.11E-02 


[CI] 


609.2 


2.75E-04 


1.77E-02 


CO 


647.2 


1.60E-03 


7.05E-01 


CO 


863.0 


1.07E-03 


5.94E-01 
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